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The overvoltage for the discharge of hydrogen on a mercury cathode has been measured at 
25°C in sulfuric acid solutions in the concentration range 0.125 to 9.65 molar. In dilute solution, 
the overvoltage is independent of the acid concentration, the cathode reaction being the 
discharge of a proton from water. At concentrations in excess of 2.83 molar, the overvoltage 
decreases approximately as a linear function of the weight percentage. Evidence is presented 
to show that in this region, the cathode reaction involves the hydrogen sulfate molecule, 
H2SQ,. Preliminary results on phosphoric acid solutions also show a drop in the overvoltage as 
the concentration is increased beyond 1 molar. 





INTRODUCTION 


HE current density of discharge of positive 

ions at a cathode is given by the theory of 

absolute reaction rates by an equation of the 
form 


ekT l= AHit 


when the process is activation controlled.! AS? 
and AH? are the entropy and heat of activation 
(per ion), respectively, for the discharge process, 
[A+] is the activity of the ions in a layer of 
adsorbed ions adjacent to the cathode, « is the 
charge on a proton, k Boltzmann’s constant, T 
the absolute temperature, and h Planck’s con- 
stant. V is the electrical potential between the 
adsorbed ion layer and its electrical image in 
the electrode; its sign is defined by the relation: 


V=V (metal) — V (solution). (2) 


* Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 
1G. E. Kimball, J]. Chem. Phys. 8, 199 (1940). 
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a is the fraction of this double-layer potential 
operating between the adsorbed ion layer and 
the top of the potential energy barrier. It is 
convenient to define the cathodic overvoltage 


n as: 
V= Vea (3) 


where V,q is the double-layer potential at an 
electrode at equilibrium in the given solution 
with respect to the given electrode reaction. 

At hydrogen cathodes in moderately acid and 
neutral solutions (0< pH <7), the experimentally 
observed current densities can be represented as 
a function of the overvoltage by an equation of 
the form: 


2.3 
I=I1,10"* =I, exp | (Tafel’s equation). (4) 


Bowden? showed that, between 0° and 72°C, the 
factor 6 does have the temperature dependence 
predicted by Eq. (1), viz.: 


. 2.3kT 59.1mv 





(at 25°C), (5) 


ae a 


2F. P. Bowden, Proc. Roy. Soc. Al26, 107 (1929). 
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with a= on most metals. This value of a was 


confirmed by numerous workers. The factor J» is 
very dependent on the metal used, depends on 
the temperature with energies of activation of 
the order of 10 to 20 kcal.mole™,’ is proportional 
to the }th power of the hydrogen pressure at 
pressures not too far below 1 atmosphere,‘ and 
is independent of the pH of the solution.® 

The theoretical interpretation of these facts 
has been given by Eyring, Glasstone, and 
Laidler,* and Kimball, Glasstone, and Glassner,? 
who showed that Eq. (1) can be made to fit the 
facts if it is assumed that (1) the source of the 
protons discharged at the cathode is water, 
(2) that the electrical double layer as it extends 
outward from the electrode into the solution 
includes at least two layers of adsorbed ions, so 
that a proton crossing the double layer must pass 
over at least two energy barriers before it reaches 
the electrode. KGG showed that the charge 
density in the inner ion layer increases with 
[H*], py,, and n, but that at reasonably high 


values of these variables, it becomes essentially 
“‘saturated,’’ and then shields the inner half of 
the double layer from any further electrostatic 
disturbances in the outer half, and vice versa. 
That part of the reversible potential which 
varies with the »H of the solution is set up 
across the outer half of the double layer, while 
that part which varies with hydrogen pressure 
operates across the inner half. If the slow step is 
the passage of a proton across the inner energy 
barrier, then the overvoltage is also set up across 
the inner half of the double layer. Thus if the 
total potential be written as the sum of an 
inner and an outer potential term :* . 


V= Vit Vo, (6) 


8 Review by F. P. Bowden and J. N. Agar, Ann. Reports 
Chem. Soc. 35, 90 (1938). 

‘S. J. Bircher and W: D. Harkins, J. Am. Chem. Soc. 
45, ~ (1923). M. Knobel, J. Am. Chem. Soc. 46, 2751 
(1924). , 

5F. P. Bowden, Trans. Faraday Soc. 24, 473 (1928). 

*H. Eyring, S. Glasstone, and K. J. Laidler, J. Chem. 
Phys. 7, 1053 (1939), hereinafter referred to as EGL. 

7G. E. Kimball, S. Glasstone, and A. Glassner, J. Chem. 
Phys. 9, 91 (1941), hereinafter referred to as KGG. 

* The standard reversible potentials V;° and V,° are 
purposely omitted from these equations, since for the 
purposes of Eq. (1), they are essentially constants which 
can be incorporated into the heat of activation AH}. 
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where 
_ ht 
Vi=— In py,—1, (7) 
2e 
and 
kT 
V,=— In (H*), (8) 


€ 


only the potential aV; is operative in activating 
the discharge reaction. From Eq. (1), the current 
density should be given by: 


—eaV; 
I=const. exp | 
kT 


= /2 | . | (9) 
const. p ailie  ' ’ 
He 





in agreement with experiment when a=}. 

The assumption that the protons come from 
water molecules accounts for the constancy of 
the factor [A+] in Eq. (1) when the fH is 
varied. The slow step, as postulated by EGL, is 
the transfer of a proton from a water molecule 
in the inner ion layer to a water molecule ad- 
sorbed on the negatively charged electrode, with 
the consequent formation of an adsorbed H 
atom on the cathode, and the liberation of an 
hydroxyl ion HO~ in the inner ion layer. Since 
the potential energy barrier exists between two 
layers of water molecules, it should be sym- 
metrical satisfying the requirement that a=3. 

The heat of activation AH? depends largely on 
the metal used. The higher the energy level of an 
H atom adsorbed on the metal (one of the 
products of the postulated slow-step), the higher 
AH}, and the greater the overvoltage for a given 
current density. Metals that are good adsorbents 
of hydrogen are known to be “‘low overvoltage”’ 
metals, and those that adsorb it poorly are 
“high overvoltage’”’ metals. By using the values 
of AH? as determined experimentally (from the 
temperature variation of I at 7 constant, or of 7 
at I constant), EGL showed that ASt is a con- 
stant independent of the metal or solution whose 
value is in fair agreement with the value calcu- 
lated for the postulated mechanism. 

However, if water ceases to be the principal 
component of the solution, as is the case when 
one goes to acid solutions of higher and higher 
concentrations, the overvoltage (J) should no 
longer be independent of the nature of the 
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solution, and should be expected to drop as 
water is progressively replaced by a stronger 
acid, i.e., by a better source of protons. This was 
actually found to be the case with both sulfuric 
and phosphoric acids, although the overvoltage 
curves (I) obtained by going to stronger acids 
were markedly different in the two cases. 


EXPERIMENTAL 


The cell used and the wiring diagram are 
sketched in Fig. 1. The cell was made of Pyrex 
glass and consisted of three compartments, for 
the anode A, the cathode C, and the reference 
electrode E, respectively. The cathode and anode 
compartments were connected by an ungreased 


l4V 
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Bridge 


Fic. 1. 


stopcock wetted with the solution being elec- 
trolyzed. This stopcock was normally kept shut 
to prevent admission of oxygen into the cathode 
compartment. The resistance of the system was 
such that currents up to about 1 ma could be 
obtained from the seven two-volt batteries in 
series (except with the most dilute solution 
employed). By opening the stopcock, currents 
up to about 10 ma could be obtained. However, 
after opening and closing the stopcock, a marked 
depression in the overvoltage 7(ZJ) amounting to 
as much as 100 mv could be observed. This was 
followed by a slow and gradual return of the 
overvoltage to its steady value before opening 
the stopcock. 








The cell was equipped with two side-arms indi- 
cated in Fig. 2. One of these led from a mercury 
funnel H through a lightly greased stopcock into 
the bottom of the cathode compartment. The 
other one led into the bottom of the reference 
electrode compartment and was used to fill the 
cell with solution. 

The reference electrode consisted of a platinized 
platinum wire sealed into a soft glass tube and 
held in place by a rubber stopper which had been 
previously boiled in acid and alkali. The platinum 
electrode was saturated with hydrogen by being 
made cathodic in situ, and then continued to 
function as a hydrogen electrode as long as the 
solution remained saturated with hydrogen. Its 
potential was controlled against that of a satu- 
rated calomel electrode K, and remained con- 
stant within better than 5 mv in the great 
majority of the runs. 

Electrical contact with the mercury cathode 
was made through a platinum wire sealed in soft 
glass, and held in place by a boiled rubber 
stopper. In this stopper was also inserted a 
stopcock which led (through a safety trap and a 
water bubbler not shown in the diagram) to the 
atmosphere and through which any excess pres- 
sure in the cell could be relieved. The area of 
the cathode was measured by running water out 
of the cathode compartment (with the platinum 
wire tube in place) and weighing the amount of 
water corresponding to a vertical distance of 1 
cm. The area was found to be 2.52+0.05 cm? 
(logio A =0.401 +0.009). 

The sulfuric acid solutions were prepared by 
dilution of 18M product (Baker and Adamson 
C.P. grade) with distilled water. Before each 
run, the cell was cleaned with concentrated nitric 
acid, sulfuric acid—K>sCr.O; solution, distilled 
water, and then steamed out. The platinum 
anode and the two platinum wire tubes were 
then inserted in place, and the stopcocks (un- 
greased) in the H»SO, train of Fig. 2 were 
moistened with some of the solution to be elec- 
trolyzed. The stopcock at the bottom of the 
cathode compartment was greased rather thor- 
oughly. The sulfuric acid solution, having been 
prepared to the approximate strength desired, 
was brought to the boiling point, and then in- 
serted in the apparatus at F (Fig. 2). The 
pressure over the sulfuric acid was reduced by 
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suction from a water aspirator, and the liquid 
allowed to boil under reduced pressure for about 
10 minutes to remove the air present. The suction 
was then cut off and hydrogen (electrolytic tank 
hydrogen, purified over red hot copper and 
through NaOH solution and soda-lime) cautiously 
allowed to bubble into the solution until the 
pressure had built up again. The hydrogen was 
then allowed to continue bubbling through the 
solution and into the electrolytic cell to displace 
the air present. When the liquid was cooled down, 
a portion of the solution was forced into the cell 
to rinse it, and discarded. A second portion of the 
liquid was then forced into the cell, hydrogen 
allowed to bubble through it for a few minutes 
longer, and then shut off. The platinum reference 
electrode was made cathodic for a few minutes 
to saturate it with hydrogen. The electrolysis 
was started between the platinum wire in the 
cathode compartment and the platinum anode. 
Redistilled mercury was quickly run into the 
bottom of the cathode compartment up to the 
level of the side arm leading to the reference 
electrode, and the overvoltage determined as a 
function of the current which could be varied by 
adjustment of the rheostat R of Fig. 1. 

In view of the depolarizing action of oxygen 
on hydrogen cathodes, the precautions recom- 
mended by all workers in this field for the 
stringent removal of oxygen cannot be too 
strongly emphasized. In the present series of 
experiments, the overvoltage »(J) generally rose 
during the first few hours of the run, gradually 
reaching a maximum and steady value. This 
indicated the presence of last traces of oxygen in 
the solution which were gradually removed by 
reduction at the cathode. This same type of 
behavior could be reproduced (as noted before) 
by opening the anode stopcock for a few minutes. 
The bubbling of ‘“‘purified’’ tank hydrogen 
through the solution led to very unsteady over- 
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voltage values. After reaching a steady value, the 
overvoltage at any given current 7(J) remained 
constant within a range generally not exceeding 
20 mv during a period of several days in which 
the measurements were taken. 

The overvoltage was taken as the potential 
difference between the platinum reference elec- 
trode and the mercury cathode: 


n= VPt(He) — Vue, (10) 


as measured by means of a Leeds and North- 
rup student potentiometer No. 152168, whose 
smallest scale division is 1 mv. The ammeters 
used were two multi-range ammeters incorpo- 
rated in radio-testing-boxes: Model 862 and 
Model 852 of the Precision Apparatus Company, 
Brooklyn, New York. In the region of interest, 
the former of these had scales of 0-12 ma and 
0-1.2 ma divided into 60 parts, and the latter 
had scales of 0-1 ma, 0-100 wa, and 0-50 ya 
divided into 50 parts. These instruments were 
used without calibration. They agreed with 
each other. 

In the majority of the runs, the apparatus 
was immersed up to line W (Fig. 1) in a water 
thermostat maintained at 25.00+0.03°C. How- 
ever, for some of the solutions, the measurements 
were repeated with the cell immersed in air at 
room temperatures of 24 to 28°C. No difference 
could be found between the two types of runs, 
since the effect of this temperature variation is 
smaller than the experimental uncertainty in the 
absolute value of the overvoltage. 

The concentration of the sulfuric acid was 
determined by titration with standard alkali on 
samples of the solution removed from the flask 
and from the cell both at the start and at the 
close of a run. The pH was determined from the 
potential of the platinum reference electrode 
with respect to calomel, and also from that of a 
glass electrode with respect to calomel (Beckman 
pH meter). The procedure used on the three 
runs with phosphoric acid reported herein differs 
in no wise from that described above for sul- 
furic acid. 


RESULTS—SULFURIC ACID SOLUTIONS 


The overvoltage measurements obtained on 


eight sulfuric acid solutions, with concentrations 
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TABLE I. Overvoltage measurements in sulfuric acid at 26°C. 








Potential Potential 








calomel to glass to log Io Range of 
Weight hydrogen calomel b (Jo in measurements 

Molarity % mv pH mv pH mv a. amp./cm?) log I 
0.125 Las 291.7+0.7 0.96 400.5 0.92 116 = 0.51 — 10.40 — 6.40 to —4.15 
0.465+0.005 44 260.1+0.3 0.42 430.7 0.41 118 0.50 — 10.30+0.12 —5.75 to —3.00 
2.35 +0.01 20.2 213.2+1.1 —0.37 116 0.51 — 10.43+0.10 —5.00 to —4.25 
130 0.45 — 9.80+0.20 —4.25 to —3.50 
2.83 +0.02 23.7 208.7+0.3 —0.45 485.9 —0.52 115 0.515 — 10.58+0.08 —6.00 to —4.25 
147 0.40 — 9.18+0.10 —4.25 to —3.00 
4.16 +0.03 32.8 191.9+0.3 —0.73 500.7 0.77 117 0.505 — 9.95+0.05 — 6.00 to —3.00 
5.85 +0.1 43.3 185.3+0.3 —0.84 $23.7 —1.16 83 0.71 — 11.70 — 6.40 to —5.25 
88 0.67 — 11.34 —5.25 to —4.50 
97 0.61 — 10.72 —4.25 to —3.50 
116 0.51 — 9.53 — 3.50 to —3.00 
8.2 +0.1 55.5 162.2+0.2 —1.23 86 0.69 — 11.06+0.03 — 6.40 to —3.00 
9.65 +0.05 62.2 138.8+0.5 — 1.63 78 0.76 — 11.50+0.08 — 6.40 to —3.00 








ranging from 0.125M (1.2 percent) to 9.65M 
(62.2 percent) are summarized in Table I and in 
Figs. 3a and 3b. If the logarithm of the current 
density is plotted against the overvoltage, straight 
lines are obtained on all solutions (except the 
5.85M) from which the constants Jp and b (or a) 
in Tafel’s equation can be determined. The con- 
stant b is given by the slope of the line and is 
equal to the rise in overvoltage accompanying a 
tenfold rise in current density; log Ip is equal to 
the intercept on the log J axis obtained by ex- 
trapolating the line to zero overvoltage. As an 
illustration of the data, the Tafel plot of the 
experimental points for runs Nos. 1 and 3 on 
the 5.85M solution is given in Fig. 4 (the other 
runs made on this solution either coincided with 
the above, or fell between the two curves and 
are omitted from the plot). This plot is not a 
straight line, but has a slight concavity towards 
the overvoltage axis. 

In Table I, the first column gives the mo- 
larity and the second the weight percentage of 
the H2SO, solutions electrolyzed. The third 
column gives the potential difference: V(calomel) 
— Vpt(H2), and the fourth the pH calculated from 
this potential (based on the value 471 mv for pH 
4.00 as measured in 1/20M potassium hydrogen 
phthalate buffer), the fifth the potential differ- 
ence: V(glass)— V(calomel) as measured in a 
Beckman pH meter and the sixth the corre- 


sponding pH value. The seventh, eighth, and 
ninth columns give the constants b, a, and log Io 
determined from the Tafel plot, and the tenth 
the range of log J within which these constants 
were observed. The values of } are obtained from 
the Tafel plot of the averaged values of the over- 
voltage, 9(J), at intervals of 0.25 in log J. The 
precision measures assigned to log Ip correspond 
to the spread in the overvoltage n{J), and are 
based on the assumption that the value of } 
chosen from the slope is free of error. Actually 
on the 2.35M and 2.83M solutions this assump- 
tion is not justified at higher current densities 
(log J>—4.25), where the Tafel plot slopes 
upwards with a comparatively large scatter in 
the values of (J), and } values ranging from 
115 to 150 mv are obtained on individual runs. 
On the 5.85M solution the Tafel plot curves 
upwards (see Fig. 4) and no single value of the 
constants can be assigned. The curve can be 
arbitrarily divided into short segments, each of 
which can be treated as a straight line. The 6 and 
log Ip values given for this solution were obtained 
in this manner. 

In Figs. 3a and 3b, the overvoltage (J) is 
plotted as a function of the sulfuric acid concen- 
tration expressed in terms of weight percent, for 
constant values of log J at intervals of 0.25 in 
log I in the range —6.25 to —3.00. The points 
represent the values of »(J) obtained from the 
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Fics. 3a and 3b. Overvoltage at constant values of the logarithm of the current 
density in sulfuric acid solutions at 25°C. 


several runs on each solution. The horizontal bar 
gives the position of the average of each cluster. 

In dilute sulfuric acid solutions, the over- 
voltage 7(J) is independent of the acid concen- 
tration and of the pH, in agreement with the 
observations of F. P. Bowden on buffer solutions 
and on a dilute sulfuric acid solution in the pH 
range 0.8 to 6.6.5 The value a=} is in agreement 
with Bowden’s, while his value of log J) on 
mercury (—11.2 at 25°C) is lower than that 
obtained here (—10.4). The overvoltage con- 
tinues to obey these characteristics as the H2SO, 
concentration is increased up to about 24 percent, 
at which point a slight peak (of ca. 15-20 mv) 
in the curve is noticeable. This peak becomes 
more pronounced at the higher current densities 
where a values below 3 are obtained, and con- 
siderable unsteadiness in the 7(J) values is ob- 
served (similar to that observed by Bowden in 
the high current density region, and attributed 
by him to concentration polarization). 

Beyond the “‘critical concentration” of ca. 24 


percent »(J) drops sharply. As the concentration 
is increased further, n(Z) continues to decrease 
approximately as a linear function of the“weight 
percent, in the region 33 percent to 62 percent. 
This drop in the overvoltage is accompanied by 
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TABLE I]. Overvoltage measurements in phosphoric acid at 25°C. 








Potential Potential 





calomel to glass to log Jo Range of 
Weight hydrogen calomel (Jo in measurements 
Molarity % mv pH mv pH b a amp./cm?) log I 

1.145+0.010 10.6 283.9+0.8 0.83 0.75 130 0.455 — 10.32+0.15 — 6.40 to —4.00 
3.21 +0.01 27.1 253.1+0.3 0.31 0.13 152 0.39 — 9.28+0.08 — 5.50 to —3.60 
5.65 +0.03 43.3 218.5+1.3 —0.28 479.3 —0.41 226 0.26 — 7.17+0.05 — 6.40 to —4,00 

191 0.31 — 7.71 —4.00 to —2.45 

266 0.22 — 6.23 —2.45 to —1.90 








a progressive rise in a from 3 to %, while log J 


remains approximately constant at ca. —11+}3. 


RESULTS—PHOSPHORIC ACID SOLUTIONS 


Overvoltage measurements were made in three 
phosphoric acid solutions with concentrations of 
1.145M (10.6 percent), 3.21M (27.1 percent), 
and 5.63M (43.3 percent). The results are sum- 
marized in Table II and Figs. 5a and 5b. The 
Tafel plot gives straight lines in all cases, and the 
constants 6 and log J» obtained therefrom in the 
manner described above for sulfuric acid, are 
given in Table II. The data on the most con- 
centrated of these solutions are represented sche- 
matically in Fig. 6 over the whole range of the 
measurements (log J from —6.40 to —1.90), the 
straight lines give the locus of the experimental 
points for the several runs on this solution. The 
slight break in the curve at an overvoltage of 
about 750 mv corresponds to the point at which 
the stopcock was opened between the anode and 
cathode compartments. 

In Figs. 5a and 5b, the overvoltage (J) is 
plotted as a function of the weight percent of 
H;PO, for constant values of log J (in the range 
—6.25 to —3.75). Again the overvoltage (J) 
is seen to decrease as the acid concentration is 
increased, but in a manner altogether different 
from that observed in H2SO,. In H3PQ,, the 
curves get further apart instead of coming to- 
gether, corresponding to a gradual decrease in a 
from a value slightly less than 3 in the most dilute 
of the three solutions, to } in the most concen- 
trated, accompanied by a rise in log Ip from 
—10.3 to —7.2. 


DISCUSSION 


The independence of the overvoltage from the 
PH or acid concentration, as observed by Bowden 


in the pH range 0.8 to 6.6, and as required by 
the EGL mechanism which assumes that the 
protons discharged come from the solvent (water), 
has been confirmed in dilute sulfuric acid solu- 
tions, and is found to hold up to a concentration 
of ca. 24 percent by weight (5.5 mole percent, 
2.83M). There is apparently a slight peak in 
the overvoltage in the vicinity of this critical 
concentration. Beyond this point the overvoltage 
n(Z) drops sharply, and decreases approximately 
as a linear function of the weight percentage of 
sulfuric acid up to ca. 62 percent (23 mole per- 
cent), the curves of constant log J tending to 
come together. It is evident that in this region 
the EGL water mechanism breaks down, and 
since the water is progressively being replaced 
by a stronger acid, it is reasonable to expect 
that the discharge of protons should require less 


TABLE III. Percentage ionization of sulfuric acid, 
from N. R. Rao’s Raman data.* 





Molarity 3.0 6. 
5 


8.7 10.8 12.8 14.9 16.9 
% Ionization’ 93.3 85. 4.1 


8 
3 7 60.3 38.2 17.2 7.1 





* There are some minor discrepancies between Rao’s reported in- 
tensities and his results calculated from them, presumably owing to 
faulty proofreading. The numbers given in Table III are taken from 
the first and last columns of Rao’s Table I. They agree with the points 
on a 2 and are, therefore, presumed to have been correctly 
reported. 


TaBLE IV. Concentration of hydrogen sulfate [H2SO,] 
and total concentration of water in the solutions used for 
overvoltage measurements. 











[H2SO,] Total concentration 

Molarity % Un-ionized moles/liter of HzO moles/liter 
0.125 0.27 0.00034 55.1 
0.465 1.005 0.00468 54.5 
2.35 5.08 0.119 50.6 
2.83 6.105 0.172 49.4 
4.16 9.00 0.374 46.4 
5.85 12.65 0.74 42.0 
8.2 22.5 1.85 35.2 
9.65 31.5 3.04 31.9 
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Fics. 5a and 5b. Overvoltage at constant values of the logarithm of the current 
density in phosphoric acid solutions at 25°C. 


activation energy, i.e., the overvoltage at a given 
current. density n(J) should drop. The prelimi- 
nary results on phosphoric acid also show that 
there is a breakdown of the water mechanism as 
the H3PQO, concentration is increased, as indi- 
cated by a drop in the overvoltage. Here, how- 
ever, the curves of constant log J tend to spread 
out further apa.t. 

It is easy to generalize Eq. (1) to the case in 
which several acids compete in donating protons 
to the cathode, and the process is activation con- 
trolled for each. It may be assumed that the 
total discharge current is given by summing 
together the rates of proton discharge from the 
different acids present, i.e., 
ekT (— AH# 


I[=— > exp | —--——_ 
i k kT 


ti aisV; HB; ’ 11 
er I + 


where the index 7 refers to the acids HB;, and 
the index j refers to the several potential terms V; 
entering into the total double-layer potential V, 
and a;; is the fraction of the potential V; that is 
operative in activating the discharge of a proton 
from the acid HB,;. Since J is given here by a 


sum of exponentials, its logarithm should no 
longer necessarily be a straight line function of 
the potential or of the overvoltage. However, as 
the potential is varied, and the different terms 
entering into Eq. (11) assume different values, 
it will generally be true that one of the terms is 
very much larger than the others in which case 
the Tafel relation still holds. 

In the system H,O-H2SOx., the possible acids 
are H.O, H;0+, HSO;-, and H.SO,. In dilute 
solutions, H,O is the main proton donor accord- 
ing to the EGL mechanism. It seems unreason- 
able to expect any appreciable amounts of the 
negative ion HSO; in the adsorbed ion layer 
adjacent to a negative electrode, and therefore 
this ion is excluded from consideration as a 
possible source of protons. The concentration of 
the ions H;O+ in the adsorbed ion layer is equal 
to the charge density in that layer, and has been 
shown by KGG to increase with [H*], py,, and 


in such a manner as to reach a “saturated”’ 
value, which accounts for the independence of 
the overvoltage (J) on that part of the re- 
versible potential drop which depends on the pH. 
Since the overvoltage does not depend on pH in 
solutions having a pH as high as 6, there is very 
strong indication that this electrically saturated 
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layer of H;O* ions is already existent at the 
cathode even in the presence of such a low hydro- 
gen ion concentration. Further increases in the 
hydrogen ion concentration in the body of the 
solution will produce only a negligible rise in 
the concentration of H,;O* in the adsorbed ion 
layer. If in these solutions of relatively high pH, 
the oxonium ions contribute only a negligible 
fraction of the discharge current J, they will not 
contribute any more to the discharge current 
when the 6H of the solution is lowered. There 
remains the acid H2SO,. The break in the over- 
voltage when the acid concentration is increased 
must, therefore, be attributed to the contribution 
of the acid H2SO, to the discharge current. The 
values of a (3 to ?) and log J) (—11+4) obtained 
in the high concentration region are then charac- 
teristic for the discharge reaction involving 
H.SO, in the presence of water as a diluent. 
This hypothesis can be tested by means of the 
Raman data on H2SO,-H.0 solutions. Of the five 
principal bands occurring in the sulfuric acid 
spectrum, it has been possible to attribute three 
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Fic. 6. Overvoltage versus the logarithm of the current 


density in 5.65M H;PQ,. The straight line segments give 
the locus of the experimental points. 


to definite molecular species (Av 910 cm to 
H2SO,, 980 to SO, and 1043 to HSO,-).8 
N. Rajeswara Rao’s’ measurements of the in- 
tensity of these three bands in solutions of 
different molarities make it possible to determine 
the actual concentration of the three species 
H.SO,, HSO,-, and SO; in a solution of any 
given molarity. His values of the “percentage 
ionization”’ (i.e., the sum of the concentrations 
of HSO,;- and SO; as a fraction of the molarity) 
are reproduced in Table III. Unfortunately, Rao 
does not give the temperature of his solutions, 
beyond stating that the cell was cooled in a 
stream of running water during the long ex- 
posures (5 hours). A smooth curve was drawn 
through the points of Table III and extended to 
100 percent ionization at zero molarity. From 
this curve, the percentage ionization was read 
off for each one of the solutions used in the 
present research, and the actual concentration of 
unionized sulfuric acid (hydrogen sulfate) calcu- 
lated (Table IV). The last column gives the 
stoichiometric concentration of water as calcu- 
lated from the weight percent and the density of 
the solution, and includes the water taken up 
in the formation of oxonium ions H,;O0*. 

If the total current density is given by the 
sum of two terms, one caused by the discharge 
of protons from H,O, the other caused by the 
discharge of protons from H2SQu,, as postulated 
by Eq. (11), then at constant overvoltage, the 
dependence of J on the composition should be 
expressible by an equation of the form: 


I=k,[H.0]+k2[H:SO,], (12) 


where the k’s are constant at a given overvoltage. 
If the water term is small compared to the 
hydrogen sulfate term, this is tantamount to 
saying that the current density is proportional 
to the concentration of hydrogen sulfate in the 
solution. That this is actually the case is shown 
by Figs. 7 to 10, which give a plot of the current 
density I as a function of [H2SO,] (from Table 
IV) for overvoltages of 509, 568, 627, and 686 mv, 
respectively. J is proportional to hydrogen sulfate 
concentration except in the dilute solutions 


8J. H. Hibben, The Raman Effect and Its Chemical 
A pplications (Reinhold Publishing Corporation, New York, 
1939), p. 356 ff. 

*N. Rajeswara Rao, Ind. J. Phys. 14, 143 (1940). 
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Fics. 7-10. Current density versus concentration of hydrogen sulfate at 
overvoltages of 509, 568, 627, and 686 millivolts. 


where the water current becomes an appreciable 
fraction of the whole, and J curves upwards 
away from the straight line. In the 4.16M solu- 
tion the currents from the two sources are 
roughly equal. 

It is therefore fairly conclusively proven that 
the discharge reaction in the more concentrated 


solutions involves the hydrogen sulfate molecule. 
This reaction has been interpreted in the above 
argument as the discharge of a proton from 
H.2SO,. The possibility that the reaction may 
actually be the reduction of HeSO, to H2SO; has 
not been eliminated either theoretically or experi- 
mentally, and may account for the markedly 
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different a and log J) values in the concentrated 
H2SO, and H3PQ, solutions. 

In the system H2O-H;PQx,, the possible acids 
are H,O, H,;0t, H.PO;,, HPO,, and H3POx,. 
The ions can be excluded from consideration for 
the reasons previously given for sulfuric acid. 
The lowering of the overvoltage in the concen- 
trated solutions can therefore be attributed to 
the discharge of protons from the acid H3PQ,. 
This process is characterized by values of a@ less 


than 4 (3 to 3) and log J» larger than those for 
water (—10 to —7). This is entirely reasonable 
since the energy of activation should be lower 
for H3PO, than for H.O corresponding to the 
higher value of log J». At the same time, if the 
potential energy curve for a proton in the H3;PO, 
molecule lies higher than that for a proton in 
HO, the activated state should move outward, 
away from the electrode corresponding to the 
lower value of a. 
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Amphipathic Character of Proteins and Certain Lyophile Colloids as 
Indicated by Absorption Spectra of Dyes* 


S. E. SHEPPARD AND A. L. GEDDES 
Research Laboratories, Eastman Kodak Company, Rochester, New York 


(Received November 18, 1944) 


A characteristic change of spectral absorption of certain dyes as between water and organic 
solvents is shown to appear in aqueous solutions of gelatin or of the colloidal detergent, cetyl 
pyridinium chloride. This behavior is attributed to a bilateral character of lamellar micelles 
of such bodies, the elementary lamellae being polar and hydrophile on one side, non-polar and 
organophile on the reverse side. The possibility is indicated of diagnosing such amphipathic 
characteristics by comparative spectrophotometric measurements with suitable dyes. 


ANY dyes, both acid and basic, and of all 
types and classes! give in aqueous solution 

an absorption spectrum differing considerably 
from that in organic solvents. The absorption in 


«exw 





TB «io* 3 x107* 


Fic. 1. Spectral absorption of Pinacyanol. Left side, in 
methyl alcohol. Right side, in water; O—4.65X 10-*M, 
x—1.0X10“*M, A—1.0X 10M, O—1.0X10-*M. 


* Communication No. 1000 from the Kodak Research 
Laboratories. 

1S. E. Sheppard, Rev. Mod. Phys. 14, 303 (1942); 
S. E. Sheppard and A. L. Geddes, J. Am. Chem. Soc. 66, 
1995-2002 (1944.) 


water does not follow Beer’s law, and the dis- 
tribution of intensity of absorption depends upon 
the concentration and the temperature. A charac- 
teristic and outstanding example is given by the 
cyanine dye, Pinacyanol, of which the alcoholic 
and the aqueous absorption spectra are shown 
in Fig. 1. 

In the course of an investigation on the diffu- 
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Fic. 2. Spectral absorption of Pinacyanol. Concentration 
=2.7X10-*M, temperature = 25°C. O—in aqueous solu- 
tion—1 hr., x—in 0.6 percent aqueous gelatin solution— 
1 hr., A—in aqueous solution—1 week, D—in 0.6 percent 
aqueous gelatin solution—1 week. 
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TABLE I. Effect of gelatin on extinctions of band 
maxima of Pinacyanol chloride solutions—1 hour old. 
Temperature, 25°C. 








«X10 
0.6 percent aqueous 
gelatin solution 


Dye concen- Aqueous solution 





tration (M) a 8B Y a B 7 

9.0 10-4 0.15 0.72 1.07 0.30 0.66 0.95 
2.7 10-4 0.22 0.87 1.04 0.32 0.67 0.73 
1.0 10-4 0.40 0.92 0.87 0.40 0.72 0.62 
1.0 10- 0.98 0.95 — 0.88 0.76 — 








TABLE II. Effect of gelatin on extinctions of band 
maxima of Pinacyanol chloride solutions—1 week old. 
Temperature, 25°C. 








«€X107 
0.6 percent aqueous 
gelatin solution 


Dye concen- Aqueous solution 


tration (M) a B ¥ a 





B ¥ 
0.9 1074 0.16 0.69 1.07 0.50 0.59 0.65 
2.4 10-4 0.23 0.75 0.97 0.51 0.63 0.47 
0.9 10-5 0.37 0.82 0.80 0.58 0.64 0.40 








TABLE III. Effect of C.P.C. on extinctions of band maxima 


of Pinacyanol chloride solutions—1 hour old. Temperature, 
75°C. 











«X10 

0.6 percent aqueous 
Dye concen- Aqueous solution C.P.C. solution 
tration (M) a 8B ¥ a B Y 
1.0 10-6 1.38 0.77 — 1.38 0.68 — 
1.0 107-> 0.95 0.95 —_ 1.06 0.66 —_ 
1.0 1074 0.40 0.93 0.88 0.59 0.69 0.52 
2.7X 10-4 0.26 0.87 1.04 0.38 0.63 0.55 
9.0 10-4 0.15 0.73 1.07 0.27 0.56 0.65 
1.8X 10-3 oo ~~ 1.10 — 0.54 0.75 








sion of dyes and electrolytes, we had occasion to 
determine the absorption spectrum of Pinacyanol 
in aqueous gelatin solutions. In the presence of 
gelatin, the absorption spectrum deviates from 
that in pure water, at the same concentration of 
dye. Furthermore, the deviation is such that the 
spectrum approaches that in organic solvents 
(cf. Fig. 1). This behavior is illustrated in Tables 
I and II, and in Fig. 2 (at one concentration 
only). 

It is apparent that the presence of gelatin 
favors the “organic’’ spectrum ¢€,>¢s>e,, and 
that there is an ‘‘aging’’ factor associated with 
this change in the presence of gelatin. (The 
change with time in the simple aqueous solution 
merely indicates a slight total fading of the dye.) 
Such a redistribution of the spectrum would be 
expected if the gelatin were so (incompletely) 


SHEPPARD AND A. 


L. GEDDES 


dispersed that its particles afforded, in part at 
any rate, a non-polar organic aspect in which 
adsorbed dye could dissociate as in a layer of 
organic solvent. The possibility of an electro- 
chemical combination of the dye and the 
(charged) gelatin? was avoided by using an acid- 
process gelatin of high isoelectric point (pH =ca. 
8.8). Consequently, at the pH of operation 
(pH=ca. 5.6), both the dye (cation) and the 
gelatin would be positively charged, and there 
would be no electrochemical (electrostatic) at- 
traction between them. Reasoning in this way, 
it appeared to us that a similar behavior to that 
of gelatin should be shown by such micelle- 
forming bodies as macromolecular detergents 
and spreading agents. Data obtained with the 
same dye (Pinacyanol) in aqueous solutions of 
cetyl pyridinium chloride (C.P.C.) appear to 
confirm the idea. See Table III and Fig. 3. 

It is evident that the cetyl pyridinium com- 
pound behaves in a manner similar to that of 
gelatin, and, for the same percentage concentra- 
tion,* rather more strongly. In both cases the 
addition of a micelle-forming hydrophile colloid 
effects conversion of some of the dye from an 
aggregated—probably dimeric—to a monomeric 
state. 

This points to the partition of the dye between 
two phases, water, and micelles of the colloid. 
At first sight it might seem that a sufficient 
picture of the process involved would be to 
assume a polar non-polar “‘amphipathic’’ dye 
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Fic. 3. Spectral absorption of Pinacyanol. Concentration 
=2.7X10-!M, temperature = 25°C. O—in aqueous solu- 
tion—1 hr., x—in aqueous solution 0.6 percent C.P.C. 
solution—1 hr., A—in aqueous solution 0.6 percent gelatin 
solution—1 hr. 


2S. E. Sheppard, R. C. Houck, and C. Dittmar, J. Phys. 


Chem. 46, 158 (1942). 
3 As yet, they cannot be compared on a molar basis. 
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molecule or ion, symbolically expressed by +40, 
and that the non-polar, organophile end pene- 
trates an oil-like globule or micelle of the 
colloid, so that “‘dimerization’”’ of the dye ions 
is prevented. However, on consideration of the 
structure of most dye molecules, it appears that 
this image is geometrically inadequate. The dye 
molecules are not linear, but planar structures; 
one edge of the planar unit may be regarded as 
containing one (or more) ionic (or polar) centers, 
and is thereby hydrophile, while the remainder 
is essentially non-polar and organophile.* Of this 
amphipathic configuration, it is the dorsal, or- 
ganophile edge or plane which would have to 
be attracted to, and cohere with, a corresponding 
portion of the colloid. 

This appears to us to be consistent with a 
lamellar rather than a spheroidal configuration 
for a colloid micelle, whether of gelatin or of 
the detergent.’ As supporting this conclusion, it 
may be pointed out that both the dye and the 


4S. E. Sheppard, Atti del X Congr. Intern. di Chimica 
(Rome) 1, 234 (1938). 
5D. G. Dervichian, J. Chem. Phys. 11, 236 (1943). 





colloid or detergent are cations. If the colloid 
particles were spheroidal, the outer surface would 
be expected to be hydrophile and positively 
charged. This would produce a strong potential 
barrier hindering the dye from penetrating the 
colloid,? whereas with a lamellar two-faced con- 
figuration, the common organophile portions can 
readily cohere. 

As more spectroscopic data for varying types 
of dye molecules (of differing constitutions, sizes, 
and shapes) become available, a new tool for 
the investigation of the micellar organization of 
colloids will be at hand. 

There are certain interesting contacts and 
relations between the phenomena just described 
and those of the ‘‘solubilization”’ of dyes recently 
discussed by McBain‘ and his collaborators, and 
also with the adsorption of dyes and proteins to 
silver halides.? It is hoped to present further 
evidence shortly on these aspects of the subject. 


6 J. W. McBain and K. E. Johnson, J. Am. Chem. Soc, 
66, 9 (1944). 

7S. E. Sheppard, R. H. Lambert, and R. D. Walker, 
J. Chem. Phys. 7, 265 (1939). 
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A Generalization of the Quasi-Chemical Method in the Statistical 
Theory of Superlattices 


C. N. Yanc* 
National Tsing Hua University, Kunming, China 


(Received November 17, 1944) 


The quasi-chemical method in the investigation of the equilibrium distribution of atoms 
in the pairs of neighboring sites in a superlattice is generalized by considering groups con- 
taining large numbers of sites. The generalized method may be used to obtain successive 
approximations of the free energy of the crystal. The labor of integration is avoided by the 
introduction of a Legendre transformation. In order to analyze the fundamental assumption 
underlying the method more closely, the number of arrangements of the atoms for given 
long-distance order is calculated and the hypothesis of the non-interference of local configura- 
tions discussed. The method is applied to the calculation of the free energy in the different 
approximations discussed in this paper, including Bethe’s second approximation and a simple 
approximation for the face-centered cubic crystal CusAu. 





1. INTRODUCTION 


T was shown by Fowler and Guggenheim! that 
the quasi-chemical method, originally devised 
for the theory of regular solutions, applies equally 
well to the theory of superlattices with long- 
distance order. The method is, as they have 
emphasized, definitely one stage further towards 
an exact theory than Bragg-Williams’ method.? 
Compared with Bethe’s* or Kirkwood’s‘ method 
it also distinguishes itself in mathematical sim- 
plicity. But in its original form the method does 
not lead to a consistent scheme of successive 
approximations and cannot be regarded as a 
method of approach towards the rigorous evalua- 
tion of the configurational free energy of an 
alloy. In the present paper a new formulation of 
the quasi-chemical method is developed which is 
capable of yielding successively higher approxi- 
mations. 

The free energy expression in both Bethe’s 
method and the quasi-chemical method involves 
an integral. Its evaluation is very complicated 
and has been carried out!* so far only in Bethe’s 


* Research Fellow of the China Foundation for the 
Promotion of Education and Culture. 

1R. H. Fowler and E. A. Guggenheim, Proc. Roy. Soc. 
A174, 189 (1940), 

2 W. L. Bragg and E. J. Williams, Proc. Roy. Soc. A145, 
699 (1934); 151, 540 (1935); 152, 231 (1935). 

3H. A. Bethe, Proc. Roy. Soc. A150, 552 (1935). 

4J. G. Kirkwood, J. Chem. Phys. 6, 70 (1938). 

5 T. S. Chang, Proc. Camb. Phil. Soc. 35, 265 (1939) ; 
J. G. Kirkwood, J. Chem. Phys. 8, 623 (1940); J. S. Wang, 
“Free energy in the statistical theory of order-disorder 
transformation,’ Science Report of National Tsing Hua 
University, Series A, 30th Anniversary Memorial Number 
(1941), printed but failed to appear. 
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approximation for simple and body-centered 
cubic crystals. In the new formulation of the 
quasi-chemical method, however, it will be shown 
that a Legendre transformation helps to effect 
the integration. (It might be noticed that a 
similar Legendre transformation is used to essen- 
tially the same effect in Fowler’s formulation of 
general statistical mechanics. Compare Fowler.®) 
The free energy is obtained directly as a closed 
expression. Its values are given for Bethe’s first 
and second (modified) approximations and for 
the face-centered alloy Cu3Au in Sections 7 and 8. 

To make sure that the quasi-chemical method 
may actually be used to obtain a series of suc- 
cessively better approximations, we must in- 
vestigate the free energy in high order approxi- 
mations and compare it with the partition 
function of the crystal. This is done in Section 5 
together with a comparison of the quasi-chemical 
and Bethe’s methods. 

Except in the last section we shall be con- 
cerned only with binary alloys with atomic 
ratio 1:1 forming a (quadratic), simple cubic, or 
body-centered cubic lattice. The generalization 
of the method to the investigation of alloys with 
other atomic ratios and forming other types of 
lattices is easy. In fact, the superior power of 
the quasi-chemical method appears to be even 
more fully revealed when a face-centered lattice 
is studied. This problem is taken up at the end 
of the paper where an approximate free energy 


6R. H. Fowler, Statistical Mechanics, second edition, 
p. 188. 
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expression for Cu;Au is obtained and its critical 
phenomena discussed. 


2. REFORMULATION OF THE QUASI- 
CHEMICAL METHOD 

Consider an alloy of type AB with altogether 
N=2N, atoms. Let z be the number of nearest 
neighbors of each atom. Those sites in the super- 
lattice occupied by the A atoms at the absolute 
zero of temperature will be called the a-sites, 
and the rest, the 6-sites. The partition function 

of the crystal may be written in the form 


> P(w, T), (1) 


where P(w, T) is equal to e~#/*? summed over 
all possible configurations of the crystal having 
Nww a-sites occupied by B (wrong) atoms. The 
average energy over all these configurations is 


_ 0 
E(w, T)=kT*— log Pw, T). (2) 
- 


N;! . 
|. (3) 
(Nw) '(N,— Nw)! 


But evidently 
P(w, ~)=g(w) -| 
Hence 


a | 
log P(w, T) =log g(w)+ f ape T)dT. (A) 





The problem of finding the partition function of 
the crystal therefore reduces to one of finding 
E(w, T). An approximate solution has been ob- 
tained by Fowler and Guggenheim. Their method 
will now be presented in a new form better suited 
to generalization. 

There are in the crystal zN,; pairs of nearest 
neighboring sites a—6. Let [ga, gg] of these be 
occupied by ga(=0, 1) wrong (B) atoms on the 
a-site and gs(=0,1) wrong (A) atoms on the 
B-site. For given w the following relations hold: 


LO, 0)+L0, 1)+([1, 0]+([1, 1]=2Mi, 
[1,0]+[1, 1]=2Niw, 


and 
[0, 1] +[1, 1]J=2Nw. 


Upon the approximation of the nearest neighbor 

interaction the energy of the crystal may be 

written as 

E(r, T)=[(0, 0) Vae+[0, 1] Vas 
+(1,0]Veet+[1,1]Vazs, (6) 


where the V’s are the interaction energies be- 
tween the different kinds of pairs of nearest 
neighbors. 

We may give (5) and (6) a different interpreta- 
tion by imagining [0,0], [0,1], [1,0], [1,1] 
and Vas, Vaa, Ves, Vaz to be, respectively, the 
numbers and the molecular internal energies of 
four different kinds of molecules, say, X, XZ, 
XY, X YZ of a gaseous assembly. The interpreta- 
tion of (6) is that the assembly has the same 
internal (non-kinetic) energy as the crystal at 
the given value of w. Equation (5) would mean 
that the total numbers of X, Y, Z atoms are, 
respectively, Ni, zNiw, 2Nyw. 

The quasi-chemical method consists in taking 
the averages (0, 0), (O, 1), (1, O)m, (1, 1) of the 
chemical assembly at any temperature as ap- 
proximately representing the corresponding aver- 
ages of the crystal at the same temperature. 

A detailed treatment of the problem of a 
gaseous assembly has been given by Fowler.’ We 
are only interested in our assembly of four 
different kinds of molecules, for which the results 
may be summarized as: 


(0, O)w=£ exp [— Vaz/kT ], 
(0, 1)w=&v exp [— Vaa/kT ], 
(1, O)w= Eu exp [—Vas/kT1, 
(1, 1)w= uv exp [— Vaz/kT ], 


(7) 


where &, u, and v are to be determined from (6). 
In Fowler and Guggenheim’s work the starting 
point is the equation 


(0, O)m(1, 1)w/(0, 1)a(1, O) av 
=exp [(Vaat Vap—2 Vap)/kT |, (8) 


which can be obtained by eliminating &, u, and v 
from our equation (7). We shall, however, make 
use of (7) instead of (8) to obtain the free energy 
by means of a Legendre transformation which 
enables us to avoid the labor of integration. It 
will appear that this procedure is applicable to 
the general case studied in this paper. 

In order to calculate the free energy from (4) 
we first write the energy of the crystal in the 


7™R. H. Fowler, Statistical Mechanics, second edition, 
pp. 162-163. 
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form 

" 0 

E(w, T) =kT’— 9, HM, Y, T), 

oT 
where 
o(é, wu, », T)=£ exp [— Vaz/kT ] 
+év exp [— Vaa/kT]+&u exp [— Vaz/kT] 
+ éuv exp [ — Vaz/kT ]. 

From (5) it is evident that 


do do dg 
&—=2M,, w—=2zNw, v—=2zNq. 
0€ Ou Ov 


If these last three equations are regarded as 
defining £, u, and »v in terms of w and T, i.e., if 
the Legendre transformation: 


log é, log uw, log », T—2Mi, 2Niw, zNiw, T 
is made, the derivative of the function 
V(w, T)=¢—2N, log E—2zNyw log p—zNww log », 
with respect to T is found to be 


OV do = 
— =— =— Ku, T). 
OT OT kT? 
This shows that the integral in (4) may be 
expressed directly as a function of &, u, vy and T: 
T 


| Be, rar=w, T)-V(w, ©). 
RT? 


to) 


The values of the parameters at T= are 
easily determined from (5) and (7): 


Ww 





é=2N,(1-—w)?, p=v= 
1—w 


Thus (4) becomes 
log P(w, T) =2 log 


N;! 
(Nw) !{Ni(i—w)}! 


2Ni(1—w)? w 
——_—+2zNww log 





+2N, log —— 
(1 —w)2yy 


so that the free energy of the crystal is 
F(w, T) = —kT log P(w, T) = —kT {2Ni log 2N, 
+2(z—1)N,[wlog w+ (1—w) log (1—w) ] 
—zN, log §¢—zNiw log wv}. (9) 


We shall have occasion to return to this ex- 
pression for the free energy later in Section 4. 


3. GENERALIZATION TO GROUPS OF FOUR SITES 


So far our attention has been fixed on the 
pairs of nearest neighbors in the crystal. They 
are classified into four different kinds and the 
average number of pairs in each class is obtained 
from chemical analogy. Now we shall generalize 
the whole procedure by taking into consideration 
all the groups of sites of a certain arbitrarily 
chosen form in the crystal. These groups will be 
classified according to the way they are occupied 
by atoms and the average number of groups in 
each class is to be obtained by chemical analogy. 

To make this clear let us consider in detail 
groups of four sites forming the corners of squares 
(as shown in Fig. 1) in a quadratic lattice. We 


Pl Rot po 
of a 


Bp 
Dapp 


Fic. 1. 


classify these groups into 24‘=16 classes and 
denote them by (0, 0, 0, 0), (0, 0, 0, 1), ---, 
(1, 1,1, 1), respectively, so that all groups in the 
class (q1, 92, 93, 94) have qi wrong atoms in their 
upper a-sites, g2 wrong atoms in their lower 
a-sites, g; wrong atoms in their upper 6-sites, and 
gs wrong atoms in their lower f-sites. The total 
number of these groups is NV, hence if we denote 
the number of groups in the class (q1, ge, gs, 44) 
by [q1, 92 da, 94], 
1 


C1, G2, 93, gs |= M1. (10) 
0 


qi= 


Now the number of all those groups in the 
crystal with a wrong atom on the upper a-site is 
just the number of wrong atoms on the a-sites. 
Similar reasoning may be applied to the other 
three sites, so we obtain 


»e giLQ1, G2, Ys, Qa |= Niw, t=1,2,3,4. (11) 


Let x(q1, 92, Ys, da) be the energy of each group 
in the class (q:, qe, qs, gs). It is easy to show that 


x= (2—qi—g2)(qs+qs) Vaa 
+ {(2—gi—qz)(2—gGs—4s) 
+ (qitqe) (93+) } Vae(2— gs—9s) (Gi +92) Vez. 


Lapa | be bo ae 
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The total energy of the crystal is given in terms 
of x(q1, q2; q3, qs) by 


E(w, T)=2 [15 G2, Ys, Ja 1x (Gis G2, Ga, G4). (12) 


We may give (10), (11), and (12) interpreta- 
tions similar to those given in Section 2 for Eqs. 
(5) and (6). The same quasi-chemical method used 
there to obtain (7) leads now to the following 
averages for given w: 


q 1y 72, Y3y 4) av = Eu Mpeg" 
Xexp [—x(q1, q2, 35 qa) /kT |. (13) 


In this expression the parameters &, 1, M2, Ms, 
and ys, are to be determined from (10) and (11), 
which may be written in the form 


dg dg 
g—=M, yw—=Nw, (1=1, 2, 3, 4) (14) 
dé Opi 
if we put 


O(E, bt, May May Ma) = u Eur po yg" 4™ 
Xexp [—x(q1, G2, 9a, Ga) /RT J. (15) 


It can be shown® that & and yu; are uniquely 
determined by (14) at given w and T. Their 
values at T= © are 


() 7-0 = Ni(i—w)!, 


w . 
(ui)r=0=——, t=1,2,3,4; (16) 
—— oe 


as can be shown by substitution into (14). 

To calculate the free energy it is necessary 
first to evaluate the integral in (4). By (12) and 
(13) the integrand may be written as 


7 1 
ae T) — x (Q1y 2» V3» G4) 


dp 
Xx(q1, 2s 35 qa) =, (17) 
oT 


where the partial differentiation is to be taken 
with £, w1, “2, us, ws, and T as the independent 
variables. If, however, we regard them as func- 
tions of w and T defined by (14) and introduce 


8 The proof follows easily from Lemma 2.42 in Fowler’s 
Statistical Mechanics, second edition. 


the Legendre transformation 
log &, log ui, log ue, log us, log ws 
—M,, Niw, Nw, Niw, Nw, (18) 
(17) reduces to 
sili T) ont T), (19) 
kT? oT oT 
where the function W is defined by 
vY¥=$—Mi log E-D Nyww log yi. (20) 


Substituting (19) into (4) we get 
log P(w, T) =log g(w) + V¥(w,T)—WV(w, ©), (21) 


so that the free energy may be written down 
directly as 


F(w, T) = —kT log P(w, T) 
= —kT[log g(w) + ¥(w, T)—V(w, ©) ]. (22) 


The equilibrium value # of w is obtained by 
minimizing F, so we have 


i) d 
0=—F(w, T)= - #7] — log g(w) 
aw dw 


0 0 
D Ow 


Ow 
But by (14) and (20) 


0 
V(w, T) = —-> Ni log Ki, 
Ow i 


and by (3) 
d w 
— log g(w) = —2N; log ——, 
dw 1-—w 
so by (16) 


Ww 





> log wi= —2 log 
i 1—w 


wD 
+{2 log wi |ran =2 log : : 


u=(—-). (24) 


i 1—w 





i.€., 








70 Cad. 


It will be shown in the next section that we 
may put V44= Vgzs, and Vag=0 without altering 
the specific heat of the crystal if 


V=3(Vaat Vaz) — Vaz 


is left unchanged. When this is done, ¢ will be 
symmetrical with respect to 1, we, w3,and us, and 
we conclude that all the yu’s are equal, because: 
(i) Eq. (14) has only one set of solutions; and 
(ii) if the conclusion is true, (14) becomes, with 
all uw; put equal to yp, 


0d 0d 
-—=M,, u— = 4 Nw, (25) 
0& Ou 


which does have a set of solutions in € and u. 
Now ¢ is given by 


b= éL1+4yx?+ (4y2x?+ 2y2x4) +4y%x?+y"], (26) 
where 
x=exp [—3(Vaat+Vaes—2Vaz)/RT }. (27) 


Introducing the degree of order s by s=1—2w 
and eliminating ~ from (25), we obtain 


(1+5)p4+ (2+4s) w3x?+ 2sx?(x?+ 2) yu? 
+(4s—2)x*u+(s—1)=0. (28) 
The free energy is given by (16) and (22): 
bit cond (1—w) log (1—w) 
2NikT 
+w log w—2w log u+4 log (1+4yx? 
+4y?x?+ Qu 2xt+4y3x2+ 44), (29) 


and the condition of equilibrium by (24): 


1—s\3 
n= (—) . (30) 
1+s 
To obtain the critical temperature, we expand 


(28) in powers of s and find, after identifying 
coefficients, 


1+6x?+x4 


2+2x? 


s+Ks*+--- 





log w= — 


which is the only real solution for log u. Next 
we expand (30); 


log w= —s—js?—--:. 


YANG 


At the critical value x, of x, these last two equa- 
tions have a multiple solution at s=0. Hence 


1 +-6x.?-+-x.‘ 
2+2x,? 


’ 


i.€., 


Xe= (4/5 —2)!=0.4858. 


4. GENERAL FORM OF THE QUASI- 
CHEMICAL METHOD 


Let us now take a group of any size and form. 
Let it have a a-sites, b B-sites, and y pairs of 
nearest neighbors. The procedures to obtain an 
approximate expression for the free energy of 
the crystal follow exactly the same line as in 
the special case considered in the last section. 
Equations (13), (14), and (16) are essentially 
unchanged : 


(q1, 92, °° w= Eur" p2x®---exp[—x/kT], (31) 
o= X Eu po - --exp [—x/kT], 


dd 
i—= Ni, S - Nw, (32) 
dg Mi 
and 
Ww 
(ui) r=0 =——. (33) 


But now the total number of pairs in the whole 
set of groups is Niy while the total number of 
pairs of nearest neighbors in the actual crystal 
is Nyz. Hence the energy expression (12) should 
be modified by multiplying the right-hand side 
with a factor z/y. Thus we have, in place of (17), 


io ; 
E=- mH (q1, d2, °° wx (G1 d2, °° ) 
vY qd 
Zz eT) 
iT — 00 
y oT 


And (22) becomes, with the help of (3), (20), 
and (33), 


z dy 

F(w, T) = ~=N:AT| tog wit (a+s-—) 
Y Zz 
X {(1—w) log (1—w) +w log w} 


—log é-—w Xu log ui} (35) 
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Differentiate (35) and we obtain 


0 
—F(w, T) 


ow 
zZ 1—w a+b—(2y/z) 
=—N,kT log | ar Mi) (—) | (36) 
Y i w 


where use has been made of (32). So the condition 
of equilibrium is given by 


w a+b—(2y/ 2) 
) (37) 


Mu (— 


1—w 


In actual calculations the following points may 
prove helpful: 

(i) The free energy is changed by a constant 
if Vas and Vgpz are both replaced by 3(Vaat Vee) 
— Vap, and Vap by 0. To prove this let 2; be the 
number of sites in the group neighboring to the 
site t. Let x be changed into x’ by the replace- 
ment. It is evident that x’—x = —yVas+(No. of 
BB pairs)(Vaa— Vee/2)+(No. of AA pairs) 
X(Vee— Vaa/2), and that 


DD @i— L Qs 


a-sites B-sites 


= (No. of BB pairs) —(No.of AA pairs). (38) 
Hence 


(G1, G2, °° * w= Eur" ps™- - -exp [—x/kT ] 


= £' 1/1 po/%- . “exp [—x’/kT], 
if we put 


t’=texp [—yVas/kT], 
us’ =u; exp [42:(Vaa— Veg) /2kT ] 


where the + sign or the — sign is to be taken 
according as the site 7 is an a-site or a B-site. 
We can now calculate the new free energy and 
verify the above statement. 

(ii) Sites that are symmetrically situated in the 
group have w’s equal irrespective of their nature 
if Vaa= Veep and Vap=0. This has already been 
shown in the last section for the special case 
considered there. With allowances for change of 
notation, the proof holds in general. Since the 
most troublesome part of the calculations is the 
elimination of the parameters, much might be 
gained by choosing a group with a large number 
of sites symmetrically situated. 


(iii) The free energy is a function of s*, so that 
(37) is always satisfied at © =} (i.e., long-distance 
order=0). The proof is simple when we have 
already made V44a=Voeep and Vag=0, so that 
an interchange of A and B atoms does not alter 
the energy: 


x(q, wo ** ‘)=x(1-q, 1—q, he -). 


For, putting 
e = Euio'*- 
and 
bi =1/ni, (39) 
we get 


Eu po®- --exp [—x/kT ] 
= £' uy! tye! 1-2. . “exp [— x/kT }. 
Thus if (32) is satisfied 


¥ (1g) ean’ ya! «exp [—x/RT] 
= u Eui%po®- --exp [—x/kT ] 


—2 ifm" po” --exp [—x/kT]=N(1—w); 


i.e., £’, wi’, we’, -** would be the solution of 
(32) with (1—w) substituted for w. Hence by 
(32) and (20) 


V(1—w, T)=Ni-M, log ¢’ 
“i N,(1 —w) log p= V(w, T) 


showing that 
F(w, T) = F(1—w, T). (40) 


(iv) The parameter for a corner site is always 
given by: 


1 
e=—{[x*s?+(1—s?)]$—sx} (41) 
1+s 


where x 1s given by (27) and s=1—2w, irrespective 
of the size of the group, of Vaa= Van, Van=O. By 
a corner site is meant a site that has only one 
nearest neighbor in the group. Let ¢e be the 
selective variable (parameter) of a corner site, 
and yw; that of its only neighbor in the group. 
If the corner site is dropped, a new group is 
obtained. We distinguish all quantities referring 
to this new group by a prime, and obtain at 
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once 
oq’ dq’ 
t’/—=N,, bi’ = Nw, 4=1, 2, ee (42) 
ag’ Ou,’ 


The sites of the primed group are numbered in 
the same way as in the unprimed group. Intro- 
ducing the variable x defined in (27) we may 
write 


b= 2. Eui™ po” - -exp [—x/kT ] 
= 2, £(1+ ex) wo®u3%- ° ‘exp [—x’/kT ] 
+a eleta)us exp C=x'/AT]. (43) 


Let these two terms be denoted by ¢o and 4), 
respectively. Since 

















dg 
¢=Ni, wi—=Nw, 
Oui 
we have $ 
go=Ni(1-—w), d1=Niw. (44) 
Now 
Odo €X 0¢1 €E 
e—= oo, e—= $1. 
Oe Ate Oe eEe+x 
Hence e—= Nww leads to 
de 
ex € 
Ni(1—w)+ Niw=Nw,. (45) 
1+ «x e+x 
or 
w e(e+x) 
= : (46) 
1—-w ite 


the solution of which is (41). Thus the two 
parameters w and y in the approximation dis- 
cussed in Section 2 are all equal to e. 
(v) The contribution to the free energy from a 
corner atom is such that, in the notation of (iv) 
y-1 1 
F(w, T) =——-F’(w, T)+—Fo(w, T), (47) 
Y sf 


where Fy(w, T) is the free energy obtained when 
y=1, i.e., the free energy (9) for the approximation 
discussed in Section 2. This is proved as follows. 
If we put 

1+ex 


’ bi=p1” =e wip”, s=2, 
1+ «x e+x 








gag” (48) 


it is evident from (43) that @ would become a 
function of £’, ui’, we”, satisfying the 





relations 
” dg dg ” do do 
— ag =M,, Mi (=1.— = Nw, 
og” O€ Ou,” Ou; 
4=1,2,---. (49) 


It is also evident that ¢ is the same function of 
ef pr”, pel’, — ¢' is of e py’, ys’’, --+, Now 
(42) has only one? set of solutions in ¢’ and u,’. 
Hence from (49) we infer that &=£, and 
ui’ = pi’. Thus 


£’ 1+ex 
f= ; ama, 


‘tthe e+x 





wi=p,’, 122. (50) 


(41) and (50) give the parameters y; in terms 
of u,’. Inserting them into (35) we obtain 


Z 2y 
F(w, T)= — NkT-| tog Nit+ (a+s-—) 
Y Zz 


X {(1—w) log (1—w) + w log w} 
—log ¢’+(1—w) log (1+ x) 


+w log (e+x) —wlog e—w dee log Xs 


y¥~1 Z 2 
a". T)— NAT: (1 --) 
ef b a 


X {(1—w) log (1—w)+w log w} 
+x 
|: (51) 


€ 





+(1-—w)log(1+ ex) +e 


If the original group (unprimed) is a pair of 
neighboring sites, we have y= 1, and (51) reduces 
to Eq. (9), i.e., to the free energy in the approxi- 
mation discussed in Section 2: 


Fo(w, T) = -svt7| (1--) 


X {(1—w) log (1—w) +w log w} 


e+x 


| (52) 





+(1—w) log (1+ €x)w log 


€ 


Inserting (52) back into (51) we get (47). 
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5. COMPARISON WITH BETHE’S METHOD 


Both the quasi-chemical method and Bethe’s 
method start from some assumption regarding 
the relative probabilities of finding a definite 
group of sites occupied in different ways by 
atoms. Although in the two methods the argu- 
ments leading to the assumptions are in no way 
similar the assumptions themselves are closely 
related. In fact, a comparison of Bethe’s local 
partition function and our function ¢ shows that 
the quasi-chemical method would give the same 
probabilities of occurrence of the local configura- 
tions as Bethe’s method if the free energy (35) 
has a minimum when 


(1) interior sites — 1 ’ (53) 


i.e., if (37) and (53) are mathematically equiva- 
lent. This is not true in general. But it holds 
approximately when the group of sites under 
consideration becomes very large. For then, 


2y 
a+b-——<y, 


Z 


hence if (53) is true, we should have 


L—wy 2th (@y/2) 717 
fay 
“ WwW 


=[]] (1) interior sites }/721 ’ 
so that by (37) 


0 
—F(w, T)=0. . 
Ow 


To see how the equilibrium free energy F(w, T) 
varies with T in high order approximations, we 
substitute (37) into (35) and obtain by making 
use of (32) 


F 1 
wn =— log (X pits”: --exp [—x/kT]) 
zNikT vy q 





1 27 
+-(a+s-—) log (1—1). 
7 oS 


“ 


The last term is very small for large groups, so 
that by (53) 


1 
=—] — x/kT }). 
war 7 og (> exp [—x/kT]) 





This shows that in high order approximations 
the free energy in the quasi-chemical method 
reduces to the exact form demanded by statistical 
mechanics. 

Recently a modification of Bethe’s method 
has been developed by Wang.’ This modified 
Bethe’s method bears a very close relation to the 
quasi-chemical method. The difference between 
the two lies in the calculation of energy which is 
discussed in length in Wang’s paper. 


6. THE NON-INTERFERENCE OF 
LOCAL CONFIGURATIONS 


Let us return to the fundamental assumption 
of the quasi-chemical method, i.e., to (31) which 
gives the average numbers of the different local 
configurations (so far called groups) in the 
crystal. This equation expresses the exact dis- 
tribution law of an assembly of molecules (cf. the 
example in Section 2) which has an energy y/z 
times as large as the crystal. Distinguishing all 
quantities referring to the assembly of molecules 
by a subscript m, we get from (4), 


1 
—F(w, T)+1 
<P, T) +log e(w) 
“| Fa T) +1 cw) 
=—-| — fp, Oo m\W . 
RT gs 


Y 


But if H is the number of arrangements in 
the crystal lattice having the given values of 


(q, lil 4 








F(w, T) = —kT log (H) +E. (54) 
Thus = = 
a. = (55) 
og ——=- log , 
(w) y gm(w 
But* . 
Ni! 
Hn = 
II qi, g2, °° - |! 


On substituting this into (55) and dropping the 
bars we get 
N,! 2/7 


H=h(w) (56) 
I Li, q2, --- ]! 





9J.S. Wang, “Approximate partition function in general- 
ized Bethe’s theory of superlattices’’ (to be published). 

*R. H. Fowler, Statistical Mechanics, second edition, 
2.6 and 5.11. 











c. 





where 


h(w) =g(w)/{gm(w) }*/7. (57) 


Equation (56) has been referred to in Fowler and 
Guggenheim’s paper! as the mathematical ex- 
pression of the ‘“‘hypothesis of the ncn-inter- 
ference of local configurations (pairs) ;’ because 
when 2/y=1, the number of arrangements in 
the crystal consistent with the distribution law 





[q1, 92, :*+ | for the groups of sites is, except for 
the factor h(r), equal to 
N;! 
H,,.= . 
at Lai, 92, --* ]! 


which would be the number of arrangements in 
the crystal having the given values of [q1, g2, --- | 
if the N,; groups in the crystal were separated and 
were filled independently with atoms. The term 
‘“‘non-interference’’ comes from the fact that the 
N, groups are not separated but are interlocked 
and cannot be filled independently with. atoms, 
i.e., they ‘interfere’ with each other. 

To find the value of g(w) we notice that by 
definition gn=2H,. But 2H, is the number of 
arrangements in the separated groups considered 
above if they are to be so filled with atoms that 
N,w of them have wrong atoms on the sites 7, 
where 7=1, 2, ---. For the N, sites 7 of the N; 
groups, 


Ni! 
[Wa(1—w)} (Wi) | 





different arrangements are possible. Hence 











N;! 
m= DL Hn= DL 
. IT La, g2, «+ J! 
N;! a+b 
-| . (58) 
Thus 
M! 2—(2/7) (a+b) 
h(w) = 59 
(w) a 1(N iw) | _ 


The free energy of the crystal may be obtained 
from (40), (42), and (45) as: 


ms 
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_ 3NkT 2y 
F(w, T)=E- | (o+0-— 
y Z 





X[(1—w) log (1—w) +w log w]+log N; 
1 f 
~z wet” G2, °° aw log (q1, do, °° nf 


which is reducible. to Eq. (35) obtained above. 


7. SPECIAL CONSIDERATIONS CONCERNING 
BETHE’S FIRST AND SECOND 
APPROXIMATIONS 


The First Approximation 


If an a-site together with its z nearest neighbors 
are taken as our group of interest, all the sites 
except the central one are corner sites. Hence 
their selective variables are all equal to the 
value of € given in (41). By successive applica- 
tions of (47) we see that the free energy is 
exactly Fo(w, T), a fact which has already been 
pointed by Fowler and Guggenheim.' The selec- 
tive variable of the central site is given by 
successive applications of (50): 


w sitex\? 
™ ( ) | 
1—w\ e+x 
The factor w/(1—w) is the selective variable of 


the central site when it alone forms the group. 
The equilibrium condition (38) becomes 





(60) 


Ae? =[0/(1—w) }*. (61) 
But by (60) and (45) 
A=[w/(1—w) }-“e*. 
Hence at equilibrium ; 
A=1. (62) 


Thus the approximation is completely equivalent 
to Bethe’s first approximation, as already men- 
tioned in Section 5. 


The Second Approximation 


Now consider the group of sites occurring in 
Bethe’s second approximation. According to Sec- 
tion 4 (iv), the selective variables for the corner 
sites in the second shell are all equal to the param- 
eter € given by (41). In Bethe’s original calcula- 
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tions, however, the selective variables for the 
corner sites and the medium sites are made equal, 
and are found to be different from e. Thus if we 
use his original method, Eq. (32) would not be 
satisfied. In other words, the probabilities of 
occurrence of wrong atoms in the corner and 
medium sites would be unequal. 

For simplicity we shall drop the corner sites 
and take as our group of interest the central 
site, the first shell sites and the medium sites 
with selective variables u, v, and X, respectively. 
The contribution by the corner sites can be 
included in the free energy by simple addition, 
as shown in Section 4. With the notations » and 
Znm Of Bethe’ we find 


O=EL (x+ux*)Palx, ¥,2), (63) 
where 


P,(x, ¥, Y=" DL uml (1+) ]" 


% (x2) (Le/21-D mm) 2) (4 2) (L221) (2-0) Con] 2) 


After the elimination of ~ and yu, (32) becomes 


In 


0 0 
ss {ast ( —u)Pd—P+uPa—P, || 
Or Or 


ZW = 





(g—2)(% x"Pn)(Q g*F) 


a 0 
{ar —u)P.r—Pi+uPr P| 
Ov 


v 


2 


n 


x _ (64) 





(d x"Pn) (do x*-*P,) 
The free energy is obtained from (35): 
NikT 





F(w, T)=—- [se—40-+4) 


s— 
X {(1-—w) log (1—w) +w log w} 
+(1—w) log (2 x"P,) 
+w log (> x-"P,) —zw log v 

; 2(z—2) 


a a log |. (65) 


8. APPLICATION TO THE CRYSTAL Cu;Au 


For the face-centered crystal Cus;Au, we may 
of course follow Peierls'® and take as our group 


10 R, Peierls, Proc. Roy. Soc. A154, 207 (1936). 








of interest a central site together with its twelve 
first shell neighbors. The free energy expression 
would then contain seven selective variables, © 
four of which may be easily eliminated. The 
resulting expression is very cumbersome and 
numerical calculations would be laborious. We 
therefore make a simpler approximation: The 
group is taken to be four nearest neighbors 
forming a tetrahedron. (See Fig. 2.) A little 


a 


rc) 
Fic. 2. 


geometrical consideration assures us that all such 
tetrahedrons contain an a-site (for gold atoms) 
and three #-sites (for copper atoms), an inter- 
esting conclusion showing that the tetrahedron 
might be regarded as a sort of ‘‘molecular’’ 
structure in a faced-centered lattice with atomic 
ratio 1:3. Our approximation may thus be 
reasonably expected to reveal the more important 
features of the order-disorder transformation in 
such alloys. 

A little reflection shows that the statement 
made in Section 4 (i) still holds, so that we may 
put Vaa= Veep and Vag=O0. Let i" and v be the 
parameters (for the wrong atoms) of the 6-sites 
and the a-sites, respectively. Let there be alto- 
gether N=4N, atoms. It is easy to see that 
there are 8N, groups in the crystal. When Nyw 
atoms on the a-sites are wrong, the equations 
determining the parameters are 


8Ni == EL? + 3x24 3x5? + x53 
+ v(x8§+3x5u+3x%u?+x%u%) |, (66a) 


do 
8N yw = v— = Ev (x®+3x9p+3x7u?+x%y5),  (66b) 


Ov 


and 


w w w 0g 
vi(~) +8y,(—) +8,(—) = p— 
3 3 3 Ou 


=3x*tul1+2xu+x4u?+ v(x+2u+xy?) ], (66c) 
where x is defined by (27). The energy of the 
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crystal is [cf. (34) ], 


. dg 
E=4kT*—-+ constant; 
oT 


(67) 


so that the free energy becomes [cf. (35) ]: 
F(w, T) = —kT[log g(w) +3(¢—8M1 log & 


—8N,w log v—8 Nw log pv) 7. |. 
But 


log g(w) = - xf 1») log (1—w) +w log w 


w 
+w log adieaiions log 1(3—w)/3} | 
and at T=, 





w w wy? 
y=—_, p= — _, ¢-8N(1—u)(1-—) , 
1-—w 3-—w 3 

Hence 
F(w, T) 
— = —9 log 3+4 log (81) 
N,ykT 


+6w log w+3(1—w) log (1—w) 
+3(3—w)-log (3—w) —4 log é 


—4w log u—4w log v. (68) 





Since ¢ and v can be very easily solved from (66), 
numerical calculations are quite simple. The 
value of the free energy is plotted in the accom- 
panying figure (Fig. 3). 

The equilibrium value of w is given by [cf. (36) 
and (37) ]: 





(1—w)(3—w) 


0= —3 log +4 log ur. (69) 


w 


This is always satisfied at w= 2.* Actual calcula- 


tion shows that the absolute minimum of the 


free energy is or is not at w=? according as 


x>0.2965 or x<0.2965. From the form of the 


*This is not evident from (69) directly. But if we 
divide the whole crystal into four simple cubic sublattices 
and introduce a w for each sublattice so that Nw is the 
number of A atoms on the jth sublattice, 7=1, 2, 3, 4, it is 
obvious that the free energy is symmetrical in the w’s. 
From this we infer that (69) is satisfied at w=3. 
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curve in Fig. 3 it is seen that the crystal has a 
critical temperature at which the long-distance 
order, and hence also the energy, is discontinuous. 
The critical temperature 7, and the latent heat Q 
are found to be given by 


1 
sitar” a Vaat Vez) — Vas], 


and 
Q=0.8824Ni[3(Vaat Vez) — Vas |. 


In terms of the total energy change from T=0 
to T=: 


Eo=3Ni3(Vaat Vee) — Vaz], 


these quantities become, with the usual notation 
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R=4Mk, 
T.=1.097Eo/R, 


and Q=0.2941Eo. 


In Bragg and Williams’ approximation, they are 
given by 


T.=2.19EF)/R, and Q=0.218E); 


and in Peierls’ approximation,’ they are given by 


' Te21.3E)/R, and Q20.36E,. 


It will be noticed that because of the lack of a 
free energy expression Peierls did not give the 
exact values of these quantities. 

In conclusion the author wishes to express his 
thanks to Professor J. S. Wang for valuable 
criticism and advice. 
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Previous theoretical treatments of the thermodynamics of high-polymer solutions have 
assumed that the entropy of mixing is independent of the magnitude of the heat of mixing. 
The present treatment considers energy interactions on a molecular scale and includes these in 
the configurational partition function. This leads to the derivation of expressions for the heat 
of mixing and the entropy of mixing which reflect the effect of solvent-solute interaction. 
One result of particular interest is an explicit expression for the temperature dependence of 


osmotic pressure. 





INTRODUCTION 


OLUTIONS of high-polymer molecules in 
ordinary organic solvents exhibit such enor- 
mous negative deviations from Raoult’s law and 
such peculiar behavior with respect to osmotic 
pressure that they hold a unique position in the 
thermodynamics of binary systems. 
In an ideal solution, the heat of mixing is 
zero and the entropy of mixing is given by: 


ASmix= —R ¥ n; in X; 


(1) 


where ; is the number of moles and X; is the 
mole fraction of the 7th component of the solu- 
tion. The partial molal entropy of mixing of the 
solvent, AS, is given by: 


AS,= —Rln Xi, 


which is obtained from the partial differentiation 
of (1) with respect to m. 

These relationships can be derived by means 
of statistical mechanics if one adopts the follow- 
ing molecular picture of a binary mixture con- 
taining N,4 molecules of the type A and Nz 
molecules of the type B. 

1. A given molecular arrangement of the mix- 
ture can be adequately described by specifying 
the occupants of all lattice sites in a quasi- 
lattice of (N4+Nz) sites. The total energy of 
the mixture can be unambiguously divided into 
two parts—one of which (independent of the 
molecular arrangement in the mixture) is simply 
equal to the energy of the separated compo- 
nents, and one of which (the “configuration”’ 
energy) depends upon the pattern of the molecu- 
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lar arrangement. Thus the partition function for 
the mixture can be factored into the partition 
functions for the two separated components and 
a configurational partition function which deter- 
mines the free energy of mixing. 

2. The configurational energy to a satisfactory 
approximation can be expressed as the sum of 
interaction energies of nearest neighbor pairs of 
site-occupants. 

3. The contribution of an A-B pair to the 
configurational energy is simply the arithmetic 
mean of the contributions of an A-A pair and a 
B-B pair, respectively; i.e., eap=(eaateps)/2. 
It follows that all arrangements of molecules in 
the lattice represent the same energy, and hence 
the configurational partition function reduces to 
a simple counting of all the possible modes of 
packing the two kinds of molecules into the 
lattice. 

4. The two molecules are of roughly the same 
size and shape, so that a given site of the lattice 
will be occupied by either exactly one A molecule 
or exactly one B molecule. 

Deviations from the laws of ideal solutions,! in 
the case of organic compounds, are generally 
attributed to failure of assumption (3) or (4) 
rather than of assumptions (1) and (2). 

If ean (€aat+enp)/2, the heat of solution will 
not be zero. Furthermore, a rigorous statistical 
mechanical treatment indicates that failure of 
assumption (3) will also change the entropy of 
mixing from the value given in Eq. (1). This 


1R. H. Fowler and E. A. Guggenheim, Statistical 
Thermodynamics (University Press, Cambridge, England, 
1939), p. 351. 








latter effect, however, is secondary in importance 
to the introduction of a non-zero heat of mixing. 
Indeed, Hildebrand and others have shown that 
the retention of the classical entropy of mixing 
even in the presence of a non-zero heat of mixing 
often represents (for mixtures of molecules of 
about the same size and shape) a very useful 
approximation. Mixtures for which this approxi- 
mation is justified have been called ‘‘regular.” 

A different type of deviation from Raoult’s law 
occurs when assumption (3) applies but assump- 
tion (4) fails. Fowler and Rushbrooke,? and 
Chang* ‘ treated the case of a solution of double 
sized molecules in a solvent of spherical mole- 
cules. In such a case, where one solute molecule 
occupies more than one lattice site of the mix- 
ture, the entropy of solution deviates from 
Eq. (1). The heat of mixing, however, is zero. 
Such solutions are designated as ‘‘athermal.”’ 

Athermal solutions of high polymers in low 
molecular weight solvents exhibit extremely 
large entropy deviations of this type. Meyer® 
discussed the problem qualitatively on the basis 
of a quasi-solid lattice, the cells of which may be 
occupied by either a solvent molecule or a seg- 
ment of the polymer molecule. Huggins® has ex- 
tended the work of Fowler and Rushbrooke in a 
somewhat quantitative manner to polymer solu- 
tions. Flory’ using Meyer’s model has arrived at 
practically an identical expression for the entropy 
of mixing. Miller* has improved upon the analyses 
of Chang by applying the Bethe approximation 
to evaluate the entropy of mixing of solvent 
molecules with “polymers” of two and three 
units similar in size to the solvent. By a con- 
sideration of the way in which the results de- 
pended upon the number of units he made an 
extrapolation to very long polymer molecules. 
The expressions for osmotic pressure, vapor 
pressure, and other thermodynamic properties 
given by these quantitative theories are in 
substantial agreement with the general behavior 
of high-polymer solutions. 

Flory and Huggins have extended their theories 

?R. H. Fowler and G. S. Rushbrooke, Trans. Faraday 
Soc. 33, 1272 (1937). 

3’ T. S. Chang, Proc. Camb. Phil. Soc. 35, 293 (1939). 

‘T. S. Chang; Proc. Roy. Soc. A169, 512 (1939). 

5K. H. Meyer, Zeits. f. physik. Chemie B44, 383 (1939). 

6° M. L. Huggins, J. Phys. Chem. 46, 151 (1942). 


7P. J. Flory, J. Chem. Phys. 10, 51 (1942). 
8A. R. Miller, Proc. Camb. Phil, Soc. 39, 54 (1942). 
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to the more general problem of a high-polymer 
solution with a non-zero heat of mixing by an 
approximation similar to that used by Hilde- 
brand® and Scatchard? in the case of low molecu- 
lar weight liquids. 

For an approximation to the entropy of mixing 
in the general case, they use the expression 
derived for the athermal case and further assume 
that the partial molal heat of solution of the 
solvent is proportional to the square of the 
solute concentration (which must be true for 
low concentrations). They suggest for the partial 
molal free energy of solution the expression : 


AF, =AH,—TAS,=kc?—TAS, (athermal). (2) 


In this present state of the thermodynamic 
theory of high-polymer solutions two aspects are 
inadequately treated. First, the entropy of 
mixing is calculated for the case where the heat 
of mixing is zero, and this value is then used for 
real solutions where there is a considerable heat 
of mixing. It is well recognized that this is not 
rigorously correct, since the presence of a heat 
of mixing must necessarily reflect a fundamental 
difference in energy among different configura- 
tions and hence a certain amount of inter- 
molecular clustering.* Such clustering will result 
in a change of the entropy of mixing from that 
calculated for the case of zero heat of mixing. 
Second, no adequate provision for the tempera- 
ture dependence of the osmotic pressure, vapor 
pressure, and other thermodynamic properties is 
included in the present theories. 

In this paper, the energy differences among 
configurations are taken into account at an 
earlier stage of the treatment, i.e., in the calcula- 
tion of the partition function itself. This permits 
the following treatment of the above-mentioned 
points. 


THE CONFIGURATIONAL PARTITION FUNCTION 


The model used here to represent a high- 
polymer solution is the simplest that is com- 


9 J. H. Hildebrand, Chem. Rev. 18, 315 (1935). 

10 G, Scatchard, Chem. Rev. 8, 321 (1931). 

*By “clustering” is meant any non-random distribution 
among configurations. This will include both deviations 
toward clusters of similar species (to be expected in cases 
of endothermic mixing) and deviations toward a regular 
“checkerboard” structure (to be expected in cases of 
exothermic mixing). 
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patible with the present conception of these 
solutions. The configurations in these solutions 
are considered to be represented by a quasi- 
crystalline lattice of m» sites, which can be 
alternatively occupied by a solvent molecule or a 
segment of a polymer molecule. This lattice 
contains m2 polymer molecules each composed of 
p segments and (%)— pm2) = n; solvent molecules. 
Our problem is to obtain the configurational 
partition function for such an assembly where 
interaction of the components results in a non- 
random distribution among configurations. The 
coordination number of the lattice is denoted 
by y. 

Provided that the internal degrees of freedom 
can be separated from the configurational degrees 
of freedom, the configurational partition func- 
tion, Q, is given by 


GM YL eFlhT=L geeFltT, (3) 


over all E, 
configurations 


where g is the degeneracy and E, is the con- 
figurational potential energy made up as follows: 


Ee= Mi1€11 + M 12€12 + Move. (4) 


Here, €11, €12, and €22 are the mutual potential 
energies of pairs of solvent, solvent-solute, and 
solute molecules, respectively, and the M’s are 
the corresponding numbers of such pairs. It is 
clear that ym;=2Mi1+My and that (y—2)pme2 
=2Mo2+Mi. In view of these linear relation- 
ships between My, Mi, and Mx, E, may be 
written simply as: 


E.=E) +My’, (5) 


where Ey’ = 4[ymieut (y — 2) pn2€22 }—i.e., the en- 
ergy of the separated components—and e’ is 
[€12— (€11 + €22)/2 ]. This can alternatively be 
written as: 


E.= Eot+ More, (6) 


where Eo, the energy of a hypothetical mixture 
in which there are no 2-2 nearest neighbors 
pairs; is defined as 


3Lymi— (y—2)pme Jerr t+ (y — 2) pres, 


and ¢€ as 


€11 + €22 — 2€12. 


The direct evaluation of 2 using (4) would be 
prohibitively complicated, but the use of Eq. (6) 
for the configurational potential energy permits 
the following solution of the problem. 

Let us return to the model, where ¥ is the coor- 
dination number of the lattice. If x polymer mole- 
cules have already been located in this lattice and 
we are considering the number of ways of placing 
the (x+1)th polymer molecule, we may distin- 
guish three types of lattice sites: (1) Those sites 
which are already occupied. (2) Those sites which 
are unoccupied, but areadjacent to occupied sites. 
(3) Sites surrounded only by unoccupied sites. 

The number of sites of the first type is ob- 
viously equal to px. We introduce the approxi- 
mation that the number of sites of type 2 is 
given by px(y—2). We thus picture each polymer 
molecule as being surrounded by an envelope of 
adjacent sites, p(y—2) in number, and assume 
that the total number of sites adjacent to 
polymer molecules is given by the sum of these 
individual envelopes. This will not be justified at 
high concentration, where the envelopes of differ- 
ent polymer molecules will interpenetrate. When 
such interpenetration becomes significant, our 
expression will overcount the number of sites 
of type 2. To the approximation that px(y—2) 
gives the number of sites of type 2, the remaining 
number of sites (type 3) will be given by 
ny— px(y—1). If there is solvent-solute inter- 
action, the availability of sites of type 2 is to be 
weighted by a factor e~*/*” which will make the 
occupation of these sites more or less probable 
according to the type of interaction present. The 
need for defining E, by Eq. (6) now becomes 
apparent, for the availability of type 2 sites 
must be characterized by a single energy param- 
eter in this model. 

The complexions* available for the (x+1)th 
polymer molecule are counted by considering 
first the addition of the terminal segment fol- 
lowed by successive segments. The number of 
available sites for this terminal segment is 


no— px(y—1) + px(y—2)e~**?, 


and the expected number of available unoccupied 
sites adjacent to the site occupied by the terminal 
* The word complexion is used to mean any geometrically 


possible configuration weighted according to how many 
type 2 sites it includes. 





segment of this polymer chain will be 


mo— px (y—1)-+px(y—2)e-ee? 
¥ . 
No 





For the third and successive segments, this 
quantity, which may be called a, will be 


Plan 1)+-px(y—2)e—*? 
/ . 


No 


(7) 





a=(y- 
The number of complexions available to a chain 
without shifting its terminal group is 


me— p(y —1)+px(y—2)ere# 
; 


No 


ar, (8) 





thus the total number of complexions for the 
chain in the lattice is 


V241=3LNo— px(y—1) + px(y—2)e-*? ] 





Y 
a?! (9) 
y—1 


x 


where the factor } compensates for the presumed 
indistinguishability of the two chain ends. To the 
approximation that y/y—1=1, Eq. (9) becomes, 
after introducing Eq. (7), 


P (=) 
_—s de 


X [no— px(y—1) + px(y—2)e-e*? ]?. 





(10) 


The configurational partition function may now 
be expressed as 
ne 
Il yz 
z=1 
Q =——e“Bolk?, 
Ne ! 





(11) 
For simplicity the factors $ and m! in Eqs. (9) 
and (11) will now be deleted. This is justified 
since the configurational partition function will 
be applied only to processes of mixing in which 
the contributions of these factors are identical 
for both the initial and final states and hence 
cancel out. 

This partition function can be expanded to 
yield a polynomial in e~*/*7: 


Q = e—Zol/kT >. pe t/kT, 
7 


(12) 
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where ¢; is the number of configurations pos- 
sessing j nearest-neighbor pairs of the type 2-2. 
Thus the process of counting the possible sites 
for each successive segment, weighting some 
sites by the factor e~*/*7 does result in the correct 
weighting of entire configurations. As given by 
Eq. (12), @ is identical to the original definition 
of Eq. (3). The expanded form (12) will not be 
used, however. 

To evaluate Ily,, the substitution of Eq. (10) 
is first made, yielding 


(: —_ -) 
No 


XII [to— px(y—1) + px(y—2)e-#/*#? J”. (13) 





By defining 





Feamy—1)-My-2e (14) 
(13) becomes 
y-—1 (p—1) ne 
( ) Il [no— ¥,}*. (15) 
No x 


Let us consider the evaluation of In II[m— Ynez], 
since it is always In Q that is required for later 
use. Hence: 


In [J [o— Yn2 |? 


=> In [mo- Ynsy= f In [mo— Yx ]*dx 


0 


=p In [my— Yx ]dx 
0 


= - (no— Ynez) —In no ] 


— Yn.[In (no— Yn2)—1]}. (16) 


This result may be used to obtain the final form 
of the configurational partition function after mo 
is replaced by (m1+ pm»). 


— Ep p ny+pn2— Yn2 


In Q= ——n, In 
n+ pne 








— 
+nz In (m1 + pn2)+(p—1)me2 In (y—1) —pne 
Dp n+ pn2— Ynz 
+(1-2) (om In ). (17) 
Y ni+pne 


The corresponding InQ for the two separated 
components (m; molecules of solvent and m2 
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molecules of perfectly arranged polymer) is given 
by —Eo’/kT. Hence the free energy of mixing 
solvent with perfectly arranged polymer is gov- 
erned by Eq. (17), with Eo replaced by (E»o— Ey’). 

To obtain corresponding expressions of the 
thermodynamic functions for the process of 
mixing ,; solvent molecules with m2 randomly 
entangled polymer molecules, one should sub- 
tract from (17) the sum of terms obtained by 
putting m; equal to zero. That is, for this latter 
process, the corresponding change in InQ is 
given by 


p n+ pn2— Ynz 
In 2’ = mecca | In 
4 mit pn 


pn p 
smi 
n+ pne ¥: 


nit pn2— Ynez pn2— Yne 
(in —In ). (18) 
n+ pne pne 


While this treatment is completely consistent 
with the fact that individual polymer segments 
may be in simultaneous contact with more than 
one other polymer segment, our partition func- 
tion gradually departs from the true value based 
on the model at concentrations where the average 
number of neighboring segments exceeds one, 
(i.e., where Ma, > pn2/2), and presumably begins 
to fail as Ma, approaches pn./2. For moderate 
values of e, this will occur at a volume fraction 
of about 1/(y7—2). 


PARTIAL MOLAL QUANTITIES 








—n» In 





The partial molal quantities that characterize 
the thermodynamic behavior of polymer solu- 
tions are obtained by differentiating the corre- 
sponding thermodynamic quantities for mixing 
which are in turn obtained from Eq. (18). 

For example, the free energy of mixing, AF nix, 
is equal to RT In 0’; that is, Eq. (18) multiplied 
by RT. When this result is differentiated with 
respect to m1, the following expression is obtained 
for the partial molal free energy of the solvent 
(free energy of dilution). 


No— Yno 





n2(p—1) 
+ 


No No 


aP.=RI|— In | (19) 
Y 


The expression for the partial molal entropy 
of mixing is obtained from (19) since — (@AF,/8T) 
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= AS. Thus, 
OY p Yn nyo— Yn 
jt: he 
OT Y*Lno— Ynez No 





-R|— In nNo— mee. aS 


No No 


|. 2 


The partial molal heat of mixing is found by 
substituting the two previous equations into 





AH, =AF,+TAS,, (21) 
and the result is 
m p aY 
AH,=+RT?*— — 
Y oT 
1 no Yno No 
x|— In + | (22) 
No no— Ynz 


The corresponding quantities for the solute 
may of course be obtained by differentiating 
first with respect to m2 or by the application of 
the Gibbs-Duhem relation. 


DISCUSSION 


From a consideration of the statistical effect 
of solvent-solute interaction on the partition 
function, the basic equations for the thermo- 
dynamics of non-regular high-polymer solutions 
have been derived. Two illustrations of the use 
of these equations will be considered here. In 
another paper we plan to present: (1) experi- 
mental confirmation of the present theory; (2) 
theoretical refinements; and (3) an extension to 
include solubility and fractionation. 

The variation of osmotic pressure with con- 
centration is of vital concern in the determina- 
tion of molecular weights of high polymers and 
the slope of the line obtained when 2/c is plotted 
against c, the concentration, is often used as an 
indication of how ‘“‘good” or “poor” is the 
solvent for the polymer. At low concentrations 
most osmotic pressure data can be represented 
by an equation of the type 


7 
—=A+Be. (23) 
c 


By the use of Eq. (19) and the thermodynamic 
relation, r= —AF',;/V;, where V;, is the partial 
molal volume of solvent, one obtains the follow- 
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ing equation for the osmotic pressure 
—RIT p n2(p—1) 
+ 


—lIn 
Vi No No 


no— Yn. 








r= | (24) 
If this expression is expanded as a power series 
in concentration, the first two terms yield : 

x RT RT Y 

i mei —e, 

Cc M 2 Vip?” p 
where V;, the molar volume of solute, replaces V; 
and p2 is the density of the polymer. If «=0, 
Y=p and this equation becomes identical with 
Flory’s and similar to Huggins’ for osmotic 
pressure of solutions in which there is no heat 
of mixing. The plot of Eq. (24) is given in Fig. 1 
for such an athermal case (A) where the arbi- 
trary assignments of molecular weight = 100,000, 
p=0.94, Vi=104 cc, and T=30°C have been 
made. In this graph, osmotic pressure is ex- 
pressed in units of cm of water. In ideal solutions, 
van't Hoff’s law is obeyed and a horizontal plot 
is obtained. In athermal polymer solutions, how- 
ever, the large entropy of mixing gives rise to a 
very steep slope for the 2/c versus c plot. This 
slope is rather insensitive to temperature (vary- 
ing directly as the absolute temperature). 

All measurements* or estimations of heats of 
mixing for polymer solutions have given only 
positive values for this quantity. From Eq. (22) 
it can be seen that this requires the net inter- 
action energy ¢ to be negative. Assigning to ea 
value of —120 calories and to y, the coordina- 
tion number, a value of 6, Eq. (24) gives, for 
30°C, the plot labeled C in Fig. 1, the slope of 
which is of the order magnitude often encountered 
in actual measurements. With this interaction 
energy, the slope of the plot has become quite 
temperature sensitive as is illustrated by plots 
B and D corresponding to 10° and 50°C, respec- 
tively. Thus it is seen that line A represents the 
osmotic pressure behavior of an ‘‘ideal’’ polymer 
solution. The further an actual plot is from this 
line, the larger is the heat of mixing and hence 
the temperature dependence of its slope. With 
increasing temperature, its slope changes toward 
that of line A. The scanty measurements of the 
temperature dependence of osmotic pressure 





* Except the system cellulose triacetate—tetrachloro- 
ethane. 
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thus far published and measurements now in 
progress in this laboratory confirm the general 
outline just described. It is now obvious that 
good temperature control is necessary for most 
measurements of osmotic pressure. 

For Fig. 1 to exhibit more quantitative sig- 
nificance it would be necessary to reduce the 
absolute values of the slopes by introducing, as 
Flory has done, a correction for the fact that 
a cell which would enclose an actual segment of 
a polymer molecule would be filled by several 
solvent molecules. 

As another illustration, one may consider the 
dependence of the partial molal heat of mixing 
on concentration. It is this quantity which has 
heretofore been represented as being proportional 
to the square of the concentration (10). However, 
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Gee and Treloar" who have made the only reli- 
able experimental determination of the quantity 
find that it varies more rapidly than the square 
of the concentration. 

Equation (22) may be transformed as follows: 


_ OV p[ no/m 1 Ne 
AH,=RT?— —| ———-+—_ In (:-v") 





oT Y No Y no 
L. no 
oY py n2\? 2 n2\°* 
nee a(2) se) 

oT YL2 No 3 No 

3 Ne ‘ 
+-1(=) tel (26) 

4 No 


1G, Gee and L. R. G. Treloar, Trans. Faraday Soc. 
38, 147 (1942). 
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by expanding and collecting terms. For the 
choice of constants: y=6, e=—100, V,=100, 
M=100,000 and 7=300°K this equation re- 
duces to 


AH, =200[c?+0.368c?+0.114c+---]. (27) 


Thus, at the low concentrations where our 
treatment is valid the square term of the series 
expansion is dominant. However, in this ex- 
ample, which is far from the extreme case, appre- 
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ciable contributions from the higher terms ap- 
pear. This treatment predicts that the AA, will 
increase with concentration somewhat more 
rapidly than C*. This is in agreement with the 
measurements of AH; by Gee and Trealoar™ in 
system rubber-benzene. 

An entirely different treatment of this problem 


by W. C. Orr has just appeared.* 


*W.C. Orr, Trans. Faraday Soc. (August, 1944). 
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The Tetrahedral X.YZ. Molecular Model 


Part II. Rotation-Vibration Energies 
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R. C. HERMAN 
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The appropriate quantum-mechanical ‘Hamiltonian is set up according to the method of 
Wilson and Howard for investigating rotation-vibration energies to second-order approxi- 
mation. The vibrational term G is derived explicitly as a function of the anharmonic constants, 
Coriolis parameters, etc. The elements are given for the secular determinant required for 
evaluation of the rotational term F including dependence on Coriolis interactions, centrifugal 
stretching, effective moments of inertia, etc. Selection rules are discussed. 


I. INTRODUCTION 


N Part I of this paper! the complete classical 
vibration problem of the tetrahedral X2YZ2 
molecular model has been presented, including 
derivation of the normal coordinates, cubic and 
quartic parts of the anharmonic potential func- 
tion, and components of vibrational angular mo- 
mentum as well as discussion of the valence-type 
potential function. In Part II the quantum- 
mechanical calculation of the rotation-vibration 
energies of the X2YZ2 model is carried out to 
second-order approximation. The quantum-me- 
chanical Hamiltonian is set up according to the 
method of Wilson and Howard? and, after con- 
tact transformation of the Hamiltonian according 
to the method discussed by Thomas,* Nielsen,‘ 
1W. H. Shaffer and R. C. Herman, J. Chem. Phys. 12, 
494 (1944). 
2E. B. Wilson, Jr., and J. B. Howard, J. Chem. Phys. 4, 
262 (1936). 


3L. H. Thomas, J. Chem. Phys. 10, 532, 538 (1942). 
*H. H. Nielsen, Phys. Rev. 60, 794 (1941). 


and Shaffer,’ a perturbation calculation of the 
allowed energies is made. 


II. CLASSICAL ROTATIONAL HAMILTONIAN 


The kinetic energy of rotation of the non- 
vibrating model has the form 


Tr=(1/2){Ie202+1eyo+lew,*}, (1) 


where the coordinate system xyz is the body- 
fixed system of principal axes defined in Section 
II of Part I and in terms of which the equilibrium 
products of inertia vanish; w,, w,, and w, are 
the components of the instantaneous angular 
velocity ; the equilibrium values of the moments 
of inertia are 


Ie2=2mb¢?+ 2n(do?+ ec") + Mee’, 
I.y= 2m(ag?+ bo?) + 2ner?+ Mc’, (2) 
I.2= 2may?+2nd,’, 


5 Shaffer, Thomas, and Nielsen, Phys. Rev. 56, 895, 


1051 (1939). 
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where the masses and geometrical dimensions are 
the ones defined in Table I of Part I. The energy 
expression (1) can be expressed in Harniltonian 
form as: 


Hr=(1/2){P2/Iez) +(Py?/Tey) +(P2/Lez)}; (3) 


where P,, P,, and P, are the components of the 
total angular momentum of rotation of the 
model. 

When the model vibrates, the moments of 
inertia have instantaneous values which depend 
on the normal coordinates, the instantaneous 
products of inertia are not zero but are functions 
of the normal coordinates and there is set up an 
internal angular momentum of vibration with 


components pz, p,, and p, defined in Eqs. (16) of 
Part I. For the vibrating rotating model the 
rotational Hamiltonian becomes: 


Hr=(1/2) {u22(P2—pz)*+uyy(Py—Py)? 
+-pse(P2—Ps)*+2p2y(P2—pz)(Py—Py) 
+ 2pye(Py—py)(P:— Pz) 
+2us2(P:—p:)(Pz—pz)}, (4) 


where P,, P,, and P, are the components of 
total angular momentum and the coefficients 
of inertia uz, etc., can be expressed to second- 
order approximation in terms of the dimension- 
less normal coordinates defined in Eq. (12) of 
Part I: 











Urz=Lee {1+ QL wi agit DL wx (anr—aintdet+fr)ar+ :--}, 
é h 

Byy=Ley MAA D wi git DL on (di? — dan tdvten)gnr+ ---}, 
i a 


Mee= Tes {1 +2 weg th wn"(en?—Corter+fean+---}, 

May= (eal) ME wf Marth wr '(enfn—andn—brda)gn?+---}, (5) 
j 

Mys= Ceyles) AX oegth wa'(dafn—bnen—Crer)ge+---}, 

Kez (esles) ND on fide 2 wn"(dren—Onfn—Crfrgne+:--}. 


The notation adopted in Part I will be used here, namely, the subscript # denotes any one of the 
nine normal coordinates; 7 assumes the values 1, 2, 3, 4, and denotes a coordinate of symmetry 
type Ai; j assumes the value 5 and denotes the coordinate of type A2; k assumes values 6, 7, and 
denotes coordinates of type B:; and / assumes values 8, 9, and denotes coordinates of type By. The 
coefficients @;, a;;, etc., occurring in Eqs. (5) are defined in the following equations in which 
t= (h/2nc); the quantities m,, are coefficients in Eqs. (7), (8), (9), and (10) of Part I giving the 
relations between the intermediate and normal coordinates; the quantities &j,", max’, and [y, are the 
Coriolis parameters defined in Eqs. (17) of Part I; and the equilibrium rotational constants are 
Bez=h/82°CI ez, Bey=h/82°cI ey, and Bez=h/8r'cl,,: 
a;=87-B,.| m(bo— Co) 11 — (Co +0) N2i— ndona; |, 
b;=817-*B.,y[ m(bo— Co) M11 — (Co +0) N2i+maons; |, (6a) 
a 
c;=87-1B,,[ mans; —ndona; |, 
a;=b;=c;=0; a,=b,=c,=0; a,=b,;=c,=0; 
Qix=2Bexl mime + pote? + ama? + 2y19M1 M2: — 32 — E97 |, 
bi= WBeyl mimi2+ pone? + wgt32+ 2p12M 12; — nei? — nied; (6b) 
CiX= 2Bee[ uss? + mana? = 57], 
a55= 2B.2[ 2n(ndo?/mac?) N55" — E65? — E757 |, 
bs5= 2Bey[|2nn58 — N35" — nos’ |, (6c) 


C55= 2Bes|_wstss? — 651? — 652" — f 53° — £54 |, 


fc 
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ann = 2Bez[ (us — p11) 6x? + (M7 — Me) M7K? + 2 (Mor — M12) M6eM7K— Exs? ], 
bik =2Bey_woten? + mmr? + 2yorMextre — ner — eo — Nes’ — Nea J; 
Crx=2BeeLuimen? +more? + 2puremermre — Fes? — Fx" ], 
Au=2Bez_ ssi? +wonor + 2ysomsitgr — E11? — E19” — E13" — Eva], 
bu=2Bey[(us— pr) Mar + (uo — we) M9? +2 (us9— 12) M8iM91— Mis’ J, 
Cu=2BesL miner t+ momen + 2pimsimg— f6r — 6777], 
d,;=e;=f;=0, 
ds=8B,z'B.y*(—ndonss), és=fs=0, di=ex=0, di=fi=0, 
fe=4Bea*Bey'[ {m(bo— 2c) +nbo} mex—2n(eo+Co) nr], 
€:=8B.y'Bes*[m(bo—Co)msi— (Co +Co) M91 J. 
The cross-product terms in Eqs. (5) and the corresponding coefficients in Eqs. (6) have been omitted 


for the sake of brevity since they do not ordinarily contribute to the energies up to second-order 
approximation. 


(6f ) 


III. THE QUANTUM-MECHANICAL HAMILTONIAN 


The appropriate quantum-mechanical Hamiltonian operator for the rotating vibrating tetrahedral 
X.2YZ» molecular model can be set up according to the method of Wilson and Howard in the form 


H=(1/2)00 X map(Pa—pa)(Ps— pe) + (he/2h*) d) wn(ubpru*pr) + U, (7) 
a 8 h 


where a and # denote the axes x, y, and 2; wag denotes one of the inertia coefficients defined in Eqs. (5) ; 
uw is the determinant of the wag; pa is the operator corresponding to the dimensionless normal coordi- 
nate gi; Pa is the operator corresponding to the component of total angular momentum; 
Pa=>-1 Xanfn is the operator corresponding to the a-component of vibrational angular momentum 
if Xa, is one of the quantities defined in Eqs. (17) of Part I; U is the potential energy function. 
The zero-, first-, and second-order parts of the Hamiltonian operator are 


9 


Hy= (hc/2) {i wal (pr?/h?) +9? ]+(2/h*) [B.2P2+BeyP,’?+B.:P2* ]} ’ (8a) 


h=1 


4 . 
H,= U,+ (he/h?) {> wi *q:[a;BezP2+5),BeyPy+c:BeeP 2? |t+wsqsdsBe2'Bey(P2Py+P,P:) 


t=1 


9 7 
+> wr ge Bey$B..*(PyP:-+P-Py+>d wn gif -Bee'Bez'(P:P2+P2P:) 


l=8 k=6 


—2(BesbzP2t+BeypyPyt+Bep-P:)}, (8b) 


9 
He= Uet (he/h?){¥ wi qn? (an? —aan td +f) BerP 2+ (bn — bin tert+dn)BeayP/ 


h=1 
+ (¢,? ‘ies Carter t+fi2)BeeP 2 |+ (Bezp2+Beypy+Bep/) } . (8c) 


The straightforward determination of the eigenvalues of H= H)+Hi+H; by means of a perturba- 
tion calculation is complicated because of the large number of variables involved and the fact that 
the tetrahedral X:YZ» model is an asymmetric rotator. The perturbation calculation is simplified by 
a contact transformation of the Hamiltonian of the type THT described by Thomas,’ Nielsen,‘ 
and Shaffer ;5 it is possible to transform from HH, into He all those terms which contribute only to 
the second-order energy E>. If the general matrix element of the transformed Hamiltonian diagonal 
in all the vibrational quantum numbers is computed with the angular coordinates and total angular 
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momentum operators treated as constants, the result yields (1) the purely vibrational contribution 
E, to the eigenvalues of H and (2) a modified rotational Hamiltonian He in which the rotational 
constants B,,, etc. are replaced by B,:, etc., effective in the respective vibrational states. The eigen- 
values of Hz must be found from a secular determinant. 

IV. THE VIBRATIONAL TERM 


The purely vibrational contribution E, to the eigenvalues of the vibration-rotation Hamiltonian 
can be written in the form, E,=hcG, where G is the vibrational term and is expressed in cm™. 
A detailed calculation shows that G can be expressed in the form 


G=GotDd wr(va td) +X [Gan(oat3)?+ s Gan‘ (Vat) (Var +3) J, (9) 
h=1 h hi <h 


where v, denotes the vibrational quantum number which is associated with the hth normal frequency 
w, and may assume positive integral values; Go is a small constant term which occurs in the second- 
order energy correction and the quantities G;,, and Gj, are the functions of the anharmonic constants, 
Coriolis parameters, etc., defined below in Eqs. (11). 

It is convenient in the definition of the G,, to introduce the following abbreviations: 


Danner = (@nton ton) (@n— wa — wns) (Wn — Wn Fon’) (Wr ton —war) |, 
Ean = (8x? — 3wn-?) / (401? — wn?) |, 
Fin = (2w1/ (401? — wn”) J, 
Crnrnrs = Lon?(wn? — wn? — on?) Danen/2], 
Kan = L(on/wn) + (wne/wn) J, 
Hoan = (wn / (402? — wn?) J. 


The quantities Go, Gas, and G,,, are given by the following relations, in which the system of sub- 
script notation is that which was described in Section II: 


Go= (1/4) © { (Gan tOant+cnn) — (3/8) (a2+02+c,7) + (1/4) (ab: +bic;+c,a;) —(d?+e?+fi) 
h 


(10) 


+(3/2)ynnnn} — (1/16) E [To Binrt Awww D ivi Biv er + AaonowD san Binw 
+4000 sv Birr +40 zon jr — 3H eeBise — 3H 5:8 jj: — 3H eB ees — 3H Burs} 
= (1/2) & | (st Eh) Beet (nist ned) Bey t (t+ Fe) Bel 
Gis= (1/4) {Gy sies— 15001Biss— D tulle, 
Gig= (1/4) L6ra5— DL wi E B45), 
Gir= (1/4) {67ecK%— Lu wi En Bees}, 
Gu= (1/4) {6yu0—- Lu w AE Bui} , 
Gi = (rare — 30 BieBavre — 3000 Bow eBos—Z (wir Biss Burirae) 


2 


2 2 
— Fi Biiv — Fe Bev Cori Biire}, 
2 2 
Gi5= {Vit53— 307 B i:B 556 — PF 5:8 5334+ Kigfi3Ber}, (11) 
2 2 2 
Gir= {size — 30 Bit Bici— FeBrectKeinnBey— DD CrrriBirn}, 
k’ 





ion 
nal 


ian 


y, 


(9) 


Icy 
nd- 


its, 


10) 


ib- 
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; . 2 2 2 
Ga= {viu—30BiBui-— FiBuitKutiBee— DL CriBiv}, 
1’ 
- 2 2 
Gje= {Viper t+K jréjeBerz— DL CigBjxr}, 
l 
2 2 
Ga={rsutKyanpBey— DX CrisB jer}, 
k 
2 2 
Gir= {yecutKetciBeze— dX Cyr yur}, 
j 
2 
Gere = {Yeeuk?— Circe Binns}, 
2 
Gw={yuvv—D CavBiv}. 
é 


¢ V. THE ROTATIONAL TERM F 


The effective rotational Hamiltonian Hp, which results from the contact transformation and 
subsequent diagonalization with respect to the vibrational coordinates, has the form: 


Hr=(he/h?){ByzP2+BoyP/7+B,-P.} —Hre, (12) 


in which Hee contains the terms associated with centrifugal stretching of the molecule and has 
the form: 


Hrc=(he/h’){r1P2 +12Py +73P2 +11(P2Py+PyP2)?+15(PyP:+P:P,y)?+16(P2Py+P:P 2) 
+7(P2P2+P7P27)+73(P2P2+P2P27)+79(PyP2+P2P,’)}, (13) 


where 
2 ‘ 2 , “ 2 ae - 
m1=B., > i (a ?/2w,") ’ T2=Bey > (b?/ 2w,*) ’ 73 Bis z, (c,/ 2w,*) ’ 
i t i 


T4=BerBey i (d?/2w;7) ; 75= BeyBez # (e?/2w/?) ; tTe=B.2Bez a (fi2/2w,") ’ 
j l = k 
T7=BerBey > (a,b;/2w,*) ’ T3= B.2Bes > # (a,c;/2w,") ’ T9=BiyBez pO (b,c;/2w,”). 


The effective rotational constants are, to second-order approximation, 


Beo= Bul +¥ an(0n+3)}, Boy= WIE Balt H)) Bes=Ball+¥ vl0s+4)} (14) 
where if Jan =[(3en2-twn?)/ (wr? —wn?)] and A =[3B i: +Biire + Bis; +Bice + Birr), 
a;=0-Ma?—au-+o-iA+2Be D Satu}, 
Bima b2— beta bd +2Bey E Sutead. (15a) 
vi= wi 62—cyutw-tecA +2Bes D Jas is}, 
j 
aj=wj df ayjt+2Beed Jinks), 
By= oF d?—dyt2Bey D Ijin}, (15b) 
yg=ost{ Cy +2Bez DL Tsitjids 


ap. =wr{ fi2?—air +2B.2 dO Justus} ° 
j 
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Be=on{—diz +2Bey DL Tuines}, (15c) 
vye=on 1 { fi2?—Crrt+2Bez Dd Tuan , 
t 


a1=wi{ —aut+2B.. > iui, 


Bi=wr{e? —but2Bey Seata. (15d) 
j 

yi=wr{er—cut2Be. dD Fetal. 
k 


The eigenvalues of Hr can be evaluated from a secular determinant which is a modified form of 
the one given by Wang® for the rigid asymmetric rotator. The determinant can conveniently be set 
up in terms of the usual rotator quantum numbers, J for total angular momentum, K for the z com- 
ponent, and M for the space-fixed Z component where J may assume the values 0, 1, 2, ---; K as- 
sumes the values 0, +1, +2, ---, +J; and M assumes the values 0, +1, +2, ---, +J. The secular 
determinant breaks up into (2J+1)-fold steps in each of which J and M are constant and K assumes 
its (2J+1) possible values. The non-vanishing elements of each step are of three general types, 
which have the following forms if f=J(J+1): 


(JKM|Hp| KM) =he{ Rit R2K?+R;K‘}, 
(JKM|Hp|JK+2M) =he{Rs+Rsl (K+2)2+K?]} {f-K(K+1)}{f—(K41)(K+2)}}, (16) 
(JKM|Hr|JK+4M) =hcRo{ f—K(K+1)}{ f—(K£1)(K+2)}{f—(K42)(K43)}}! 


X {f-—(K+3)(K+4)}}, 
where 


Ri=(f/4)[2(Bozt+Boy) + rit r2+4174—2( 15+ 16+ 77) + (f°/8)[—371- 372-474-277], 

R2= (1/8)[(8B,.—4(B.2+ Bry) —571—572—20(74— 75— 76) +1077] 
+(f/4)[3714+372+474—875—876+27;—473—4 179 ], 

R3= (1/8) —371—372+873—1274+16(75+ 76) —277+8(738+ 75) |, 

Ra= (1/4) [(Boz— Boy) +5(16— 75) —f(71— 72) J, 

Rs=(1/8)[71— 2 —476+475—273+275 |, 

Re=(1/16)[ — 11 — 72 +474 +277 ]. 

The eigenvalues of Hr obtained from the secular determinant can be expressed in the form Er=hcF, 

where F is the rotational term. A method has been discussed by Shaffer and Nielsen’? and Nielsen*® 


whereby it is possible to factor each JM step of the secular determinant into four steps and thereby 
facilitate the evaluation of the rotational term, F. 


VI. SELECTION RULES 


The vibrational selection rules were stated in Section II of Part I. The rotational selection rules 
are those given by Dennison? for the asymmetric rotator. Dennison has also discussed the effect of 
the identical particles on intensities of the rotational lines in the rotation-vibration bands. 


6S. C. Wang, Phys. Rev. 34, 243 (1929). 

7 W. H. Shaffer and H. H. Nielsen, Phys. Rev. 56, 188 (1939). 
8H. H. Nielsen, Phys. Rev. 59, 565 (1941). 

9D. M. Dennison, Rev. Mod. Phys. 3, 280 (1931). 
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Further Studies of the Structure and Stability of Burner Flames'~* 


GUENTHER VON ELBE AND Morris MENTSER 
Central Experiment Station, Bureau of Mines, Pittsburgh, Pennsylvania 


(Received September 23, 1944) 


It is shown that, contrary to Mache’s model of the progressive extinction of a flame from 
the burner rim, a combustion zone cannot vanish within a combustible stream. The depth of 
penetration of the quenching of the explosive reaction by the burner wall is calculated from 
values of burning velocity and critical boundary velocity gradient for flashback; it is compared 
with the limiting distance between plane-parallel] plates and the limiting tube diameter for 
flame propagation. The thermal expansion of the gas normal and parallel to the combustion 
zone is discussed. An experimental analysis and discussion of partial entrance of the combustion 
zone into the burner tube (tilted flame) and partial attachment to the burner rim are given. 
New data have been obtained on hydrogen and acetylene flames. For instantaneous flashback, 
the boundary velocity gradients are independent of tube diameter, as expected; these gradients 
are not a satisfactory criterion for flame stability because flashback can be readily induced by 
tilted flames. The limit of the tilted-flame range is represented by the semitheoretical equation 


g/(1—4S,,/gR)*=const. 


g is the boundary-velocity gradient, S, the burning velocity, and R the tube diameter. The 
boundary-velocity gradients for blow-off are again found constant over the laminar flow 
range. The compression of acetylene-oxygen streams by the combustion zone has been 
measured. The burning velocities calculated from these and additional thermodynamic data 


agree well with those determined from gas flow and cone surface. 





N a previous paper,‘ the flashback and blow-off 
of burner flames, the pattern of the gas flow, 
and the distribution of gas velocity and tempera- 
ture have been investigated theoretically and 
experimentally. The experiments were made with 
mixtures of natural gas and air. The present 
paper contains additional observations on flames 
of mixtures of hydrogen and acetylene with air 
and oxygen, and a discussion of several points 
that have not been fully treated before. 

In this, as in the previous paper, the term 
“combustion zone” is synonymous to the fre- 
quently used terms ‘‘explosion wave”’ or ‘‘flame 
front.”’ It denotes the narrow preheating and 
reaction zone which separates the burned from 
the unburned gas. Within the zone, in the layers 
adjacent to the completely burned gas, heat is 
released and atoms and free radicals are formed 
by the chemical reaction; a part of the heat and 
active species (chain carriers) is transferred by 

1 Published by permission of the Director, Bureau of 
Mines, U. S. Department of the Interior. 

2From a thesis submitted by Morris Mentser to the 
faculty of the Graduate School of the University of Pitts- 

burgh in partial fulfillment of the requirements for the 
degree of Master of Science. 

§ Contribution from the Central Experiment Station, 
Bureau of Mines, U. S. Department of the Interior. 


‘B. Lewis and G. von Elbe, J. Chem. Phys. 11, 75 
(1943). 


conduction and diffusion to the adjoining layers 
of unburned gas which in turn become activated 
to rapid chemical reaction. In this way, the 
combustion zone progresses against the unburned 
gas at a definite velocity S,, the burning velocity, 
whose magnitude is evidently affected by any 
factors which modify the reaction rate and the 
rates of heat transfer and diffusion into the un- 
burned gas. Thus, at a certain distance from a 
solid wall, the burning velocity decreases to 
zero, and this distance becomes smaller as the 
temperature of the wall is raised. 

The combustion zone is generally plainly 
visible owing to chemiluminescence, notably by 
emission of CC and CH bands; and if the zone 
retains a fixed position against a gas stream 
issuing from a burner port, it forms the well- 
known inner flame cone. 

I. SUPPLEMENTAL DISCUSSION OF THE 
STABILITY AND STRUCTURE OF 
BURNER FLAMES 
(A) The Propagation and Structure of the 
Combustion Zone 


While the work of Mache‘ has materially con- 
tributed to the theory of burner flames, it appears 


5H. Mache, Die Physik der Verbrennungserscheinungen 
(Veit und Comp., Leipzig, 1918), pp. 25-54. 
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Fic. 1. Illustration of Mache’s concept of progressive 
extinction of the combustion zone from the rim of the 
burner port. 


that neither this author nor subsequent investi- 
gators have taken sufficient account of the 
mechanism of propagation of the combustion 
zone. Thus, an erroneous concept has persisted 
in the literature regarding a supposed tendency 
of the burner flame to extinguish itself pro- 
gressively from the rim toward the stream axis. 
This is illustrated in Fig. 1. For simplicity, the 
gas velocity U, is taken as uniform over the 
cross section of the stream issuing from the 
burner port. An inner cone is supposed to have 
been formed with an opening angle 2a@ so that 
S,=U, sina; that is, the burning velocity is 
everywhere equal to the normal component of 
the gas velocity, and a necessary condition for 
the stability of the cone is thus fulfilled. It is 
now further argued that any element of the 
combustion zone (as, for example, the element at 
point O) is simultaneously carried upward by 
the gas flow with the velocity U,, and inward by 
the combustion with the velocity S, and that 
consequently the element moves along the re- 
sultant of the vectors U, and S, within the 
cone surface at a velocity U, cosa. On this 
basis it would appear, as illustrated according 
to Mache in Figs. 1b and 1c, that an element at 
the rim moves within the combustion zone 
toward the cone tip; that is, the cone is pro- 
gressively extinguished from the rim toward the 
tip. To maintain a stationary cone, Mache 
assumes that the combustion zone is continually 
restored at the cone base by a virtual pilot 
flame which is maintained by a supposed slow- 
moving, turbulent boundary layer surrounding 
the gas stream. 

To visualize the consequences of this concept, 
consider an elementary tube extending normal 
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to the combustion zone from the element O into 
the unburned gas. This tube encloses all the gas 
elements that enter the combustion zone suc- 
cessively at points above O, and the tube moves 
correspondingly upward. Since, in the absence 
of a pilot flame at the cone base, the gas at 
lower points is supposed to remain unignited, 
it is virtually assumed that the reaction-inducing 
flow of heat and chain-carriers from the suc- 
cessively burning gas elements propagates only 
within the tube. This is certainly not the case. 
Rather, the progress of a combustion zone may 
be visualized in accordance with the principles 
of Huygen’s construction for the propagation of 
a wave front, as shown schematically in Fig. 2A. 
A gas element entering the combustion zone is 
activated to rapid chemical reaction by contribu- 
tions of heat and chain-carriers, not only from 
the preceding gas element in the same stream 
tube, but also from gas elements in the adjoining 
stream tubes; in turn it activates, not only the 
succeeding gas element in its own stream tube, 
but also the elements in adjoining stream tubes. 
Thus, ignition spreads from any point of the 
combustion zone laterally as well as forward; 
accordingly, a combustion zone cannot vanish 
within a combustible stream as Fig. 1 would 
indicate; rather, it always tends to occupy the 
whole cross section of the stream and moves as a 
unit with or against the stream or remains 
stationary, as the case may be. The combustion 
zone moves with the stream—the flame blows 
off—when the gas velocity everywhere exceeds 
the burning velocity ; it moves against the stream 
—the flame flashes back—when the burning 


























Zone of heat flow | emia, 
<a emadicceen 
1 3 4 A 
Stream tube “ \ 
i 
1 | zp 
B 
Tig 
2 ihe 3 
| 











Fic. 2. A. Schematic view of the combustion zone. 
B. Temperature and concentration gradients within the 
combustion zone. 1. Concentration of fuel gas. 2. Tem- 
perature. 3. Concentration of atoms and free radicals. 
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(a) Velocity gradients at wall of burner tube. (b) Velocity gradients at boundary of stream 
Curve 1, flash-back; curve 2, flash-back above orifice. Curve 1, stable flame; 
limit; curve 3, no flash-back. curve 2, blow-off limit; curve 3, blow-off. 


Fic. 3. Illustration of flashback and blow-off conditions. 


velocity exceeds the gas velocity at any one 
point; and the flame remains stationary when 
the burning velocity equals the gas velocity at 
any one point but does not exceed it anywhere. 

Other details of the combustion zone that are 
shown schematically in Fig. 2 are self-explana- 
tory. In accordance with previous practice, the 
subscripts u and 6 refer to unburned and burned 
gas, respectively. The temperature at the be- 
ginning of the zone of active chemical reaction 
may be referred to as the ignition temperature 
Tig; it obviously bears little relation to the 
ignition temperatures determined by heating or 
adiabatically compressing a combustible mixture 
to the point of explosion because the ignition 
process is controlled by other factors besides 
temperature. 


(B) Depth of Penetration of the Quenching of 
the Explosive Reaction by the Burner Wall 


It has been shown previously that a burner 
flame remains stable between two critical values 
of the gradient of gas velocity at the stream 
boundary.‘ For laminar flow in cylindrical tubes 
the gradient g may be written as a function of 
the total gas flow V and the radius R, 


—dU 
lim ( )=2=4V/aR (1) 
r 


rR 


d 


If subscripts F and B refer to the flashback and 
blow-off condition, respectively, then gr and gz 
are the lower and upper critical values of the 
velocity gradient between which the flame is 
stable. The conditions near the stream boundary 
are illustrated in Figs. 3a and 3b. The values of 
gr and gz obviously represent the slopes of the 
gas velocity curves 2 in the corresponding dia- 
grams. With reference to Fig. 3a, the distance X 
from the wall over which the burning velocity is 


below its normal value may be defined as the 
depth of penetration of the quenching of the 
explosive reaction by the burner wall. It can be 
calculated approximately from measured values 
of gr and S, by assuming that, at the point of 
contact of the curves of burning velocity and 
critical gas velocity, the burning velocity is close 
to its normal measured value. It then follows that 


X = S,/gr. (2) 


The value X may be considered a measure for 
the limiting distance between plane parallel 
plates and, correspondingly, for the limiting 
diameter of cylindrical tubes below which flame 
propagation is quenched. The validity is only 
approximate, since X refers to a single wall 
bordering on a large volume of gas, whereas the 
quenching effects of two opposing walls merge 
and augment each other. Hence, the limiting 
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Fic. 4. Depth of penetration of the quenching of the 
explosive reaction by the burner wall. Natural gas-air 
mixtures. Curve 1, limiting diameter of cylindrical tubes 
for downward flame propagation (O data by Lewis and 
von Elbe; 0’ new data). Curve 2, limiting distance between 
plane parallel plates for downward flame propagation. 
Curve 3, depth of penetration calculated from critical 
velocity gradients for flashback. 
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Fic. 5. Diagram of vertical center plane of natural gas-air flame on rectangular 
burner tube of 0.755X2.19 cm cross section. Mixture composition, 7.5 percent natural 
gas; flow, 204 cubic centimeters per second. (a) and (b): Combustion zone and flow 


pattern. (c) Burning velocity. 


distance between plane parallel plates is larger 
than X and may, in the absence of a quantita- 
tive theory, be estimated to be comparable to 2X. 
The limiting diameter of cylindrical tubes is even 
larger because of the additional effect of curva- 
ture. This is illustrated in Fig. 4. Curve 1 
represents previously obtained data‘* with two 
new points added in order to insure that the 
composition of the natural gas had not changed 
appreciably in the interval between previous and 
present experiments. Curve 2 represents observa- 
tions on a burner formed by two plane parallel 
plates, one fixed and the other movable. Curve 3 
represents values of 2X calculated from pre- 
viously published data of gr and S,.‘ It is seen 
that the limiting distance between plane parallel 
plates is even larger than 2X, although not 
greatly so. The data of Fig. 4 refer to walls at 
room temperature; as the temperature is in- 
creased, the depth of penetration decreases. 

It is of interest to report an observation by 
F. A. Smith® that a mixture of 2H2+Oy, flashed 
back in a cold glass tube of 0.013-inch (0.033-cm) 
diameter. For this mixture gr is found to be 
1.3105 sec.—! (see Fig. 8) and S,, according to 
Jahn,’ is 870 cm per second. The value of X is, 
therefore, 0.0067 cm which accords well with 
Smith’s observation. 


® Combustion (American Gas Association, New York, 
1932), 3rd edition, p. 110. 

7G, Jahn, Der Ziindvorgang in Gasgemischen (Oldenbourg, 
Berlin, 1934), p. 21. 


(C) Modification of the Gas Flow at Tip and 
Base of the Conical Combustion Zone 


In Fig. 5, a previous analysis‘ of a flame cone 
and flow pattern is reproduced with an additional 
curve showing the change of burning velocity 
over the cone surface. The analysis shows that 
the gas elements that enter the combustion zone 
in the range of constant burning velocity between 
tip and base expand only in the direction normal 
to the zone. At the tip and base, however, 
thermal expansion also takes place parallel to 
the combustion zone. In the axial stream tube, 
the expansion appears to be almost completely 
of the latter type, since the experimental records 
of this and other flames show that there is no 
appreciable change of gas velocity between un- 
burned and burned gas in the central flow line. 
This lateral expansion at the cone tip is made 
possible by the divergence of the flow lines 
caused by refraction in the combustion zone. 
At the base the expanding gas elements acquire 
considerable velocity components parallel to the 
combustion zone in the direction of the free 
atmosphere. This is clearly shown by-the outer- 
most flow lines of Fig. 5b. Since the acceleration 
of the gas elements results in opposing inertia 
forces, the expansion parallel to the combustion 
zone in the inner flow lines is gradually checked, 
and thermal expansion normal to the zone 
predominates. 

In the original photographs, the flow lines are 
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Fic. 6. Flame of natural gas and air, partially drawn into burner tube at gas flow 
just exceeding critical flow for flashback. Mixture composition, 8.1 percent natural gas in 
air; gas flow, 52 cc per second; inside diameter of tube, 1.068 cm. A. Photograph with 
particle tracks in vertical center plane (cf. Lewis and von Elbe, reference 4). B. Diagram 
of flow pattern and velocity distribution. Solidly-drawn flow lines are so located that the 
mass flow between any two adjoining lines is the same. 


seen to diverge slightly even before entering the 
heating zone. It may be remembered that the 
combustion zone exerts a pressure on the un- 
burned gas because of the acceleration of the gas 
flow. This acceleration is smaller at the base of 
the cone than elsewhere because the burning 
velocity there is lower. Consequently, there 
exists a pressure drop in the unburned gas below 
the combustion zone from the stream center to 
the burner rim which appears to be the cause of 
the observed divergence of the flow lines. 


(D) The Tilted Flame 


A flame form, frequently observed when the 
gas flow is gradually reduced toward the flash- 
back condition, is typified by the diagram in 
Fig. 6. It has already been described briefly.‘ 
The dynamic symmetry of the gas stream is 
destroyed: The combustion zone has partly 
entered the burner tube and exerts a pressure 
normal to the tube axis which causes a deflection 


of the stream of unburned gas in the opposite 
direction. The resulting change of velocity dis- 
tribution in the unburned gas causes a decrease 
of the boundary velocity gradient, below the 
critical flashback value, in the region where the 
flame penetrates into the tube; this permits 
the combustion zone to propagate against the 
stream. As the combustion zone travels down- 
ward, the resistance of the flow of unburned gas 
to the deflecting pressure increases until a posi- 
tion of equilibrium is reached, in which the 
boundary velocity gradient at the deepest point 
of penetration has become equal to the critical 
gradient for flashback, but exceeds the critical 
value at lower points in the burner tube. Ob- 
viously, the equilibrium position is maintained 
only when the burner tube is cooled sufficiently. 
Generally, as the tube temperature increases, 
the combustion zone continues to travel down- 
ward, and flashback results. 

When the gas flow is gradually increased 
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Fic. 7. Critical velocity gradients at boundary of gas stream for flashback of 
hydrogen-air flames in cylindrical water-cooled tubes of different diameters. 


toward the blow-off condition, the combustion 
zone often remains attached to the port rim at a 
single point and extends from there at a steep 
angle over the gas stream. The phenomenon is 
analogous to the foregoing : The pressure exerted 
by the combustion zone deflects the flow lines 
in the unburned gas so that the boundary velocity 
gradient at the point of attachment is smaller 
than the critical gradient for blow-off. At all 
other points along the port rim, the critical 
gradient is exceeded. On further, usually slight, 
increase of the gas flow, the flame blows off 
completely. 

Tilted flames, such as are shown in Fig. 6, are 
formed only below some critical flow which is a 
function of tube diameter. It is obviously de- 
sirable to determine this function theoretically, 
but the problem appears to be rather difficult 
and has not been completely solved. Neverthe- 
less, it has been found possible to arrive at a 


semi-empirical equation with the aid of theoretical 
considerations. The equation is given in the 
following experimental part of the paper. 


II. EXPERIMENTS ON HYDROGEN AND 
ACETYLENE FLAMES 


(A) Apparatus and Procedure 


The two gases, comprising the gas mixture to 
be studied, were passed individually through 
calibrated flowmeters into a mixing chamber. 
Effective mixing was accomplished by directing 
the two jets in the chamber at a right angle 
against each other (see, for example, Rummel’). 
On leaving the chamber, the gas mixture entered 
directly into a vertically mounted burner tube 
whose tip was surrounded by a water jacket. 
The gas was ignited by approaching the stream 
near the rim of the tube with a small flame. 

8K. Rummel, Der Einfluss des Mischvorganges auf die 


Verbrennung von Gas und Luft in Feuerungen (Verlag 
Stahleisen, Diisseldorf, 1937). 
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For all experiments hydrogen, oxygen, and 
acetylene were supplied from tanks of com- 
pressed gas, while natural gas and air were 
taken from the lines in the laboratory. No 
purification of these gases was considered neces- 
sary, other than the removal of water vapor from 
the air by passing it through a calcium chloride 
trap. 

Flowmeters for small rates of gas flow (rates 
not exceeding 4 to 5 cc per sec.) were calibrated 
by the method of displacement of water from a 
graduated vessel. With all flowmeters of larger 
capacity a wet test meter was employed. Because 
of the comparatively high solubility of acetylene 
in water, neither of the foregoing methods of 
calibration was applicable to this gas. Therefore, 
a system was devised by which a volume of 
acetylene could be withdrawn from the tank, 
measured, and then displaced under constant 
pressure through the flowmeter by water that 
was maintained saturated with acetylene gas. 

In all calibrations the pressure at the exit of 
the flowmeter was atmospheric. During the ex- 
periments the exit pressure generally rose a few 
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Fic. 8. Critical velocity gradients at boundary of gas 
stream for flashback of hydrogen-oxygen flames in cylin- 
drical water-cooled tubes of different diameters. 
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Fic. 9. Critical velocity gradients at boundary{Yof gas 
stream for flashback of acetylene-oxygen flames in cylin- 
drical water-cooled tubes of différent diameters. 


mm Hg above atmospheric pressure because of 
the flow resistance of the mixing chamber and 
the burner tube, and corrections became neces- 
sary to reduce the flowmeter readings to volumes 
delivered at the burner port. Theoretically, if the 
flow through the meter follows Poiseuille’s law, 
the volume passing through the capillary per 
unit time at a given scale reading is independent 
of pressure; hence, the volume delivered at the 
port is increased by the factor (1+AP/P.), 
where AP is the excess of pressure above atmos- 
pheric at the flowmeter exit and P, is the. 
prevailing atmospheric pressure. If the constric- 
tion of the flowmeter consists of a plate orifice, 
the volume passing through the orifice per unit 
time at a given scale reading is proportional to 
the inverse of the square root of the pressure at 
the meter entrance (see, for example, Whitwell’) 
and hence, the volume delivered at the port is in- 
creased approximately by the factor (1+AP/Pat)! 
~1+4AP/P.. Tests have shown that, for the 
square-cut capillaries used in this work, the 
correction factor could be written as 1+yAP/Pat, 
where for air, oxygen, and acetylene, y=0.5, 
and for hydrogen, y=0.8. 

Pyrex-burner tubes were used for the series of 
experiments with hydrogen-air mixtures. For 
mixtures of hydrogen and oxygen, tubes of 


® J. C. Whitwell, Ind. Eng. Chem. 30, 1157 (1938). 
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Pyrex and a porcelain thermocouple tube with 
an inside diameter of 0.058, cm were used. The 
burner tubes for acetylene-oxygen flames were 
made by drilling holes of 0.049 to 0.071 cm 
diameter in brass rods. To provide for laminar 
gas flow at the orifice, the length of any burner 
tube was at least 40 times its inner diameter. 

Determinations of limits of inflammability be- 
tween plane parallel plates for natural gas-air 
mixtures of 6.0 to 11.5 percent natural gas were 
made with the use of a rectangular metal burner 
whose cross section was of 1.6-cm length and of 
‘adjustable width. The limiting width for the 
downward propagation of a given mixture was 
measured with a micrometer. Inflammability 
limits were redetermined in a cylindrical tube of 
known diameter and compared with older data 
to preclude errors arising from possible changes 
in natural gas composition. 

The pressures developed in the acetylene- 
oxygen flames were measured by a water barom- 
eter, which, at one side, was attached to the 
burner tube just below the orifice, and open to 
the atmosphere at the other side. A continuous 
flow of water through the jacket surrounding the 
burner tip was required in order to obtain con- 
stant values. 


(B) Ranges of Flashback and Flame Tilt 


From the values of the flow rates at which 
instantaneous—or true—flashback of hydrogen- 
air and hydrogen-oxygen flames occurred, the 
velocity gradients at the wall may be calculated. 
Plots of the critical velocity gradient versus 
mixture composition for these two types of 
flames (solid curves, Figs. 7 and 8) show the 
behavior expected theoretically, that is, the 
velocity gradient for flashback of a flame of a 
given mixture composition is independent of the 
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Fic. 10. Idealized combustion zone of tilted flame. 


TABLE I. Calculation of K at the limit of the tilted 
flame range for various gas mixtures. 








Rcm g sec.71 K X10-, sec.-! 





9 percent natural gas+91 percent air; S.=38 cm/sec. 
0.775 652 


0.78 
0.709 504 0.67 
0.649 493 0.68 
0.534 494 0.65 
0.440 433 0.96 

34 percent hydrogen+66 percent air; S,=230 cm/sec. 
0.226 18,100 21 
0.147 14,600 19 
0.108 13,500 22 
60 percent hydrogen +40 percent oxygen; S,=790 cm/sec. 

0.0622 265,000 300 
0.0381 158,000 230 
0.0292 122,000 360 








tube diameter. However, the foregoing flash- 
back data lose their practical significance be- 
cause the regions of tilted flame, shown by the 
areas between the dotted and solid curves of 
Figs. 7 and 8, become much more extensive than 
is the case with natural gas-air flames. Actually 
these regions constitute a continuation of the 
range of unstable flame, inasmuch as the develop- 
ment of a tilted flame is followed invariably by 
flashback except for very lean or rich mixtures. 

Flashback induced by flame tilt, as distin- 
guished from the instantaneous phenomenon, is 
rather easily recognized with hydrogen flames. 
While hydrogen-air flames are not distinctly 
visible because of their low luminosity, the 
existence of tilt, though short-lived, can be 
observed in a shadow picture of the flame; with 
hydrogen-oxygen flames, the brilliance and pin- 
point smallness of the cone often obscures the 
partial entrance of the flame into the tube, but a 
time-lag is observed between the formation of a 
flame and flashback. Rich flames of acetylene- 
oxygen mixtures may be drawn into the tube at 
several places around the rim, thus giving rise 
to a pronged flame which sometimes rotates 
because the points of entry of the flame shift 
from time to time. Since, as discussed below, 
the process of ignition induces flashback in 
acetylene-oxygen mixtures at flows considerably 
exceeding the true flashback limit, these mixtures 
were first ignited at flows far within the stable 
flame range, and the flow was then reduced until 
flashback occurred. In this process the flame 
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' became strongly tilted and tended to flash back flow lines, that is, the pressure that the com- 
— on slight disturbances; therefore, the flashback bustion zone exerts normal to the tube axis. 
ell limits shown in Fig. 9 probably represent the If the combustion zone is assumed to be a plane 
Cc. lowest obtainable limits induced by tilted flames, surface within the burner tube, as shown in 
and the region of true flashback cannot be Fig. 10, the component of pressure normal to 
outlined exactly. the tube axis may be taken proportional to the 


On the whole, the partial entrance of the com- sine of the angle ¢; and, at the limit of the 
bustion zone into the burner tube is a more tilted flame range, the ratio g/sing may be 
_ complicated phenomenon in acetylene-oxygen expected to remain the same for a given mixture 
flames than in other flames. The limits of the composition in tubes of different diameters. 
tilted flame range could not be reproducibly From geometrical relations 
established. They seemed to depend largely on RN} Rv} Rv} 
fortuitous changes in the contour of the burner sin ¢= (:-—) = (1+=) (1-=) , (3) 
port. The intense heat and the radial gas flow 0” o o 


—_ at the flame base (cf. Fig. 5) melted away the , , arse — 
_ ( g. 5) if where o is the major semi-axis of the elliptical 


original sharp edge of the tube and formed an : : . 
8 ial combustion zone of area moR. Since the ratio 


sec, 
























































sh- irregular recess around the port. : . : 
e- Since tilted flames are formed by partial flash- ibs —. names to 1, Eq. (9) 
he back, it appears that, above some limiting value — y 
of of the gradient g of gas velocity at the tube wall, : aaa 
an tilted flames cannot be formed. However, this sin s-14(1-—) . (4) 
ly limiting value changes with tube diameter, since 
he the phenomenon depends on the deflection of the Equating the gas flow V with the gas flow through 
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Fic. 11. Critical velocity gradients at boundary of gas stream for blow-off of hydrogen- 
air flames in cylindrical water-cooled tubes of different diameters. Arrows with symbols 
mark the beginning of turbulent flow (Reynold’s number 2200) in the tubes corresponding 
to the symbols. 






















VELOCITY GRADIENT AT BOUNDARY, # =4¥., SECONDS” ‘ 
Fic. 12. Critical velocity gradients at boundary of gas 


stream for blow-off of acetylene-oxygen flames in cylin- 
drical water-cooled tubes of different diameters. 


the combustion zone 
V=70RS,, (5) 


where 5S, is taken to be constant over the zone. 
Substituting Eq. (1) in (5) and transforming 


R/« = 4S,,/gR, (6) 


and the condition that, at the limit of the tilted 
flame range, g/sin ¢ remains constant, may be 
written 


. = 
(1—4S,,/gR)! 





(7) 


Table I contains numerical values of K for the 
various gas mixtures investigated. The data for 
natural gas-air mixtures are taken from Lewis 
and von Elbe.*?®° The burning velocities for 
hydrogen mixtures are taken from Jahn.’ The 
constancy of the K values is satisfactory. 


(C) Range of Blow-off 


Figures 11 and 12 show the critical velocity 
gradients for blow-off of hydrogen-air and acetyl- 
ene-oxygen flames. The hydrogen-oxygen flames 
have not been investigated because their blow-off 
curves generally fall outside the laminar flow 
range. This is also largely true for hydrogen-air 
flames. Since Eq. (1) for the velocity gradient is 
valid only for laminar flow, the extension of this 
equation to the range of turbulent flow yields 
values of g of obscure physical significance that 
form curves branching off the main laminar flow 
curve. By combining Eq. (1) with the equation 


10 The apparent anomalous extrusion of the gr curve for 
a tube of 1.550-cm diameter (cf. reference 4, Figs. 8a 
and 11) has been recognized to belong to the tilted-flame 
range. 
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for the Reynolds number, 
2p V/xnR= 2200, (8) 


the gradient g; at the point of transition from 
laminar to turbulent flow is found to be 


£1=2n(2200)/R?p. (9) 
It is seen that each curve of Fig. 11 undergoes a 
pronounced change in slope at the value of g, 
(indicated by arrows in the figure), computed 
for the corresponding tube. 

Figure 11 also shows that the critical velocity 
gradients for blow-off become dependent on the 
tube diameter when the diameter and the gas 
flow are very small. Under these conditions the 
quenching effect of the port rim penetrates into 
the interior of the stream beyond the linear 
range of the gas velocity curve; therefore, blow- 
off is attained with a smaller velocity gradient. 


TABLE II. Pressure difference between unburned gas and 
atmosphere for acetylene-oxygen flames. 








Acetylene in 
mixture Total gas flow Flame pressure 
percent cc/sec. cm 





Inside diameter of burner tube, 0.049 cm 


18.0 6.0 2.0 
28.0 6.0 a 
46.0 6.0 3.4 
Inside diameter of burner tube, 0.061 cm 

14.0 4.0 1.2 
8.0 0.9 
12.0 1.2 

Av. 1.1 
18.0 6.0 2.4 
9.0 2.3 
12.0 235 

Av. 2.4 
23.0 6.0 3.4 
9.0 3.4 
12.0 3.6 

Av. 3.5 
28.0 8.0 5.0 
10.0 5.3 
11.0 5.4 
12.0 Sz 
14.0 5.7 

Av. 5.3 
36.0 6.0 5.6 
9.0 5.4 
12.0 55 

Av. 5.5 
46.0 6.0 3.4 
8.0 3.3 
10.0 3.7 

Av. 3.5 
52.0 4.0 1.3 
5.0 1.5 
6.0 1.5 
8.0 1.7 

Av. 1.5 
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STRUCTURE OF BURNER FLAMES 


TABLE IIT. Calculated temperatures and density ratios for acetylene-oxygen flames of various mixture compositions. 
(T,=291°K, P=1 atmos.) 


99 



















CsH2 in mixture, 
percent co COz 2 





Mole fraction of reaction products at equilibrium To —_ 
Oo Oo OH H:0 


He H (°K) Py 








d 0.183 0.471 0.075 0.098 0.074 0.005 0.011 2980 10.75 
28.6 0.326 0.118 0.099 0.113 0.144 0.078 0.026 0.090 3303 14.61 
35.7 0.441 0.076 0.028 0.076 0.118 0.073 0.057 0.127 3366 15.95 
45.5 0.556 0.021 — 0.021 0.035 0.035 0.135 0.198 3357 18.13 
50.0 0.600 — — — — 0.206 0.194 3288 18.83 











An equivalent statement is that, for a given 
velocity gradient, flames of richer mixtures can 
be blown off in small tubes. This is illustrated by 
the lower part of Fig. 11, where the critical 
gradients for blow-off in small tubes lie above 
those determined in larger tubes. 


(D) Flame Pressure and Burning Velocity 


The pressure difference on both sides of the 
combustion zone that is caused by the mass 
acceleration within the zone is given by the 
equation® (cf. Lewis and von Elbe) "4 


Pu— Po=Ap=pwS.?(pu/ps—1). (10) 


The pressure p, represents essentially the pres- 
sure throughout the whole stream of unburned 
gas because the pressure gradient along the 
burner tube is very small; ~, is substantially 
identical with the pressure of the surrounding 
atmosphere; so that Ap may be taken as the 
pressure difference between the gas in the burner 
tube and the outer atmosphere. 

For natural gas-air mixtures the flame pressure 
Ap, calculated from values of S, and p./p»,* is 
only 0.01 cm of water for the strongest mixture ; 
on the other hand, for acetylene-oxygen mixtures, 
pressures up to 5.5 cm of water are observed. 
This is shown in Table II. The resulting com- 
pression of an acetylene-oxygen mixture in the 
burner and connecting tubes at the instant of 
ignition slows down the gas stream momentarily 
to such an extent that flashback occurs over a 
considerable range of gas flows that would 
otherwise yield stable flames. In determining the 
flashback point, it is therefore necessary to ignite 
the mixture at a much higher flow, and to reduce 
the flow gradually until flashback occurs. 





11 B. Lewis and G. von Elbe, Combustion, Flames, and Ex- 
— of Gases (Cambridge University Press, 1938), pp. 


Since the pressures of acetylene-oxygen flames 
can be measured conveniently and rather accu- 
rately, they may be used to calculate the burning 
velocities. For this purpose, p./p» is calculated 
from the equation : 

=m, T> 


Pu/ p= —, (11) 
ts Te 





where =m is the total number of moles and T the 
absolute temperature. No correction is required 
for the small pressure difference between burned 
and unburned gas. Values of 2m,/2m, and T,, 


TABLE IV. Inner cone dimensions and burning velocities 
calculated from cone areas of acetylene-oxygen flames. 








Acetylene Total Cone Cone Cone Burning 
in mixture gas flow base, height, area velocity (Su) 
percent cc/sec. mm mm mm? cc/sec. 





Inside diameter of burner tube, 0.049 cm 


18.0 6.0 0.858 1.330 1.89 318 
8.0 0.779 1.759 2.21 363 

10.0 0.771 2.083 2.57 389 

28.0 6.0 0.761 1.112 1.41 427 
8.0 0.731 1.374 1.63 490 

10.0 0.717 1.697 1.96 515 

36.0 6.0 0.782 1.149 1.49 402 
8.0 0.873 1.461 2.09 383 

10.0 0.784 1.737 2.19 456 


Inside diameter of burner tube, 0.071 cm . 


14.0 12.0 0.959 1.954 3.03 396 
16.0 0.985 2.550 4.02 398 


18.0 9.0 0.859 1.288 1.83 491 


12.0 0.945 1.588 2.46 488 
15.0 0.905 2.006 2.92 514 
23.0 12.0 1.015 1.362 2.32 518 
15.0 1.022 1.771 2.96 507 
18.0 1.078 2.090 3.70 493 
36.0 15.0 1.041 1.650 2.83 530 
20.0 1.028 2.112 3.51 570 
46.0 12.0 1.192 1.650 3.29 365 
15.0 1.116 1.999 3.64 384 
52.0 10.0 1.080 1.815 3.21 312 
13.0 1.131 2.210 4.05 321 
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Fic. 13. Ratio of densities, p, and ps, of unburned and 
burned gas and absolute temperature of acetylene-oxygen 
flames. Initial temperature 291°K, pressure =1 atmos. 


for T,=291°K and atmospheric pressure, are 
calculated thermodynamically, as described else- 
where,” with all possible dissociation equilibria 
considered. The results are summarized in Table 
III and Fig. 13. For mixtures of 50 percent to 
28.6 percent acetylene, the calculations are given 
by Montagne." Since these calculations do not 
extend to leaner mixtures, the calculation for 
the 14 percent mixture is added. 

The usual method of determining burning 
velocities, as originally proposed by Gouy, con- 
sists in measuring the area F of the inner cone 
and the gas flow V; taking the gas flow normal 
to the combustion zone as S, X F, one obtains 


S.= V/F, (12) 
where S, represents an average value from tip 


12 Reference 11, pp. 285-93. 
BM. P. Montagne, Chaleur et Saduntste 19, 54 (1938). 
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BURNING VELOCITY, CENTIMETERS PER SECOND 
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Fic. 14. Burning velocity curve for acetylene- 
oxygen mixtures. 


to base. It is hardly worth while to determine 
the cone area precisely; instead, F may be 
calculated from the cone height h and the radius 

’ of the base by the equation for the right 
circular cone 


=7R’(R?+h?)}. (13) 


Table IV contains values of 2R’ and h, measured 
by means of a cathetometer, for various acety]l- 
ene-oxygen flames. It may be noted that in 
flames of this type the diameter of the cone base 
considerably exceeds the tube diameter. Constant 
values are obtained with a burner tube of 0.071- 
cm diameter, while in the smaller tube of 0.049- 
cm diameter, the gas flow appreciably affects the 
results, indicating a considerable influence of 
the port rim. Figure 14 shows the values of S, 
obtained from flame pressures and cone dimen- 
sions, respectively; although the two methods 
are profoundly different, the results agree very 
well with each other. 





